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order moments and Renyi functions, which are nonlinear, hence displaying the multifractality 
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^SJ , 1. Introduction 

(N 

Multifractal models have been used in many applications in hydrodynamic turbulence, finance, ge- 

■ nomics, computer network traffic, etc. (see, for example, Kolmogorov f94f , 1962, Kahane 1985, 1987, 
Novikov 1994, Frisch 1995, Mandelbrot 1997, Falconer 1997, Schertzer et al. 1997, Hartc 2001, Reidi 
2003). There are many ways to construct random multifractal models ranging from simple binomial 
cascades to measures generated by branching processes and the compound Poisson process (Kahane 
1985, 1987, Falconer 1997, Schmitt and Marsan 2001, Harte 2001, Barral and Mandelbrot 2002, 2010, 

k> ; Bacry and Muzy 2003, Riedi 2003, Morters and Shich 2002, 2004, Shieh and Taylor 2002, Schmitt 

■ 2003, Schertzer et al. 1997, Barall et al. 2009, Ludena 2008, 2009). Jaffard (1999) showed that Levy 
processes (except Brownian motion and Poisson processes) are multifractal; but since the increments 
of a Levy process arc independent, this class excludes the effects of dependence structures. Moreover, 
Levy processes have a linear singularity spectrum while real data often exhibit a strictly concave 
spectrum. 

Anh, Leonenko and Shieh (2008, 2009a, b, 2010) considered multifractal products of stochastic 
processes as defined in Kahane (1985, 1987) and Mannersalo, Norros and Riedi (2002). Especially 
Anh, Leonenko and Shieh (2008) constructed multifractal processes based on products of geometric 
Ornstein-Uhlenbeck (OU) processes driven by Levy motion with inverse Gaussian or normal inverse 
Gaussian distribution. They also described the behavior of the q-th order moments and Renyi func- 
tions, which are nonlinear, hence displaying the multifractality of the processes as constructed. In 
these papers a number of scenarios were obtained for q G Q H [1, 2], where Q is a set of parameters 
of marginal distribution of an OU processes driven by Levy motion. The simulations show that for 
q outside this range, the scenarios still hold (see Anh, Leonenko, Shieh and Taufer (2010)). In this 
paper we present a rigorous proof of these results and also construct new scenarios which generalize 
those corresponding to the inverse Gaussian and normal inverse Gaussian distributions obtained in 
Anh and Leonenko (2008), Anh, Leonenko and Shieh (2008). We use the theory of OU processes with 
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tempered stable law and normal tempered stable law for their marginal distributions (see Barndorff- 
Nielsen and Shephard 2002, Terdik and Woyczynski 2004, and the references therein). Note that the 
log-tempered stable distribution is also known (up to constants) as the Vershik-Yor subordinator (see 
Donati-Martin and Yor 2006, and the references therein). 

The next section recaptures some basic results on multifractal products of stochastic processes as 
developed in Kahanc (1985, 1987) and Mannersalo, Norros and Riedi (2002). Section 3 contains the 
general Cq bounds for cumulative process of multifractal products of stationary processes. Section 
4 establishes the general results on the scaling moments of multifractal products of geometric OU 
processes in terms of the marginal distributions of OU processes and their Levy measures. Similar 
results for the finite and infinite superpositions of OU processes are proved in the Section 7. The 
number of multifractal scenarios with exact forms of the Renyi function are given in the Sections 8-13. 

Our exposition extends results of Mannersalo, Norros and Riedi (2002) on the basic properties of 
multifractal products of stochastic processes. We should also note some related results by Barndorff- 
Nielsen and Schmiegel (2004) who introduced some Levy-based spatiotemporal models for parametric 
modelling of turbulence. Log-infinitely divisible scenarios related to independently scattered random 
measures were investigated in Schmitt and Marsan (2001), Schmitt (2003), Bacry and Muzy (2003), 
Rhodes and Vargas (2010), see also their references. 



2. Multifractal products of stochastic processes 



This section recaptures some basic results on multifractal products of stochastic processes as developed 
in Kahane (1985, 1987) and Mannersalo, Norros and Riedi (2002). We provide an interpretation of their 
conditions based on the moment generating functions, which is useful for our exposition. Throughout 
the text the notation C, c is used for the generic constants which do not necessarily coincide. 
We introduce the following conditions: 

A'. Let A(i), t E R+ = [0,oo), be a measurable, separable, strictly stationary, positive stochastic 
process with EA{t) = 1. 

We call this process the mother process and consider the following setting: 
A". Let A{t) = A'^'^ i — 0,1,... be independent copies of the mother process A, and A['^ be the 
rescalcd version of A^*) : 

A['^(t) = A«(t6''), t€R+, i = 0,1,2,..., 

where the scaling parameter 6 > 1, and = denotes equality in finite-dimensional distributions. 
Moreover, in the examples, the stationary mother process satisfies the following conditions: 
A'". Let A{t) — cx-p{X{t)}, t G M+, where X (t) is a strictly stationary process, such that there 
exist a marginal probability density function 7r(a;) and a bivariate probability density function 
p{xi, X2',ti — t2). Moreover, we assume that the moment generating function 

A/(C) = Eexp{CX(t)} (1) 

and the bivariate moment generating function 

M(Ci,C2;ii -i2) = Eexp{CiX(ti) + C2X(t2)} 

exist. 

The conditions A'-A'" yield 

EAi^)(t) = A/(l) = l; 
Y&T:k^\t) = M{2) -l = al<oo- 

Cov{A['\h),A\;\t2))^M{l,l;{h-t2W)-l, b>l. 
We define the finite product processes 

n ( ^ 'I 

An(i) = n^i'^W=c^P E^^''(^^') Me [O'l]' (2) 

i=0 I i=0 J 
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and the cumulative processes 

Anit)= f A„(s)ds, n = 0,l,2,...,t G [0,1], 
Jo 

where = 0, ... ,n, are independent copies of a stationary process X{t),t > 0. 

We also consider the corresponding positive random measures defined on Borel sets B of K.-|_ 

^ln{B)= / K[s)ds, n = 0,1,2,... 

J B 



(3) 



(4) 



Kahane (1987) proved that the sequence of random measures /z„ converges weakly almost surely to 
a random measure /i. Moreover, given a finite or countable family of Borel sets Bj on ]R_|_, it holds 
that lim„_>.oo fJ-n{Bj) = fJ-{Bj) for all j with probability one. The almost sure convergence of A„ (t) in 
countably many points of R+ can be extended to all points in if the limit process A (t) is almost 
surely continuous. In this case, lim„_>.oo An(t) = A{t) with probability one for all t e R_|_. As noted in 
Kahane (1987), there are two extreme cases: (i) An{t) A{t) in Ci for each given t, in which case 
A{t) is not almost surely zero and and is said to be fully active (non-degenerate) on R+; (ii) A„(l) 
converges to almost surely, in which case A{t) is said to be degenerate on R+. Sufficient conditions 
for non-degeneracy and degeneracy in a general situation and relevant examples are provided in 
Kahane (1987) (Eqs. (18) and (19) respectively.) The condition for complete degeneracy is detailed in 
Theorem 3 of Kahane (1987). 

The Renyi function of a random measure fx, also known as the deterministic partition function, is 
defined for t G [0, 1] as 



T{q) 



lim inf 



log /, 



2"-l 



(5) 



fc=0 



where /, 
b. 



in) 



[fc2-", (k + 1)2-"] , fc = 0, 1, . . . , 2" - 1, 



An) 



is its length, and log^ is log to the base 



Remark 1. The multifractal formalism for random cascades and other multifractal processes can be 
stated in terms of the Legendrc transform of the Renyi function: 

T*(q) =min (qa — T(q)) . 
In fact, let / (a) be the Hausdorff dimension of the set 

iog/.(4")(o) _ 



Cq. 



t e [0, 1] : lim 

n^oG 



log 



where (t) is a sequence of intervals I^"^ that contain t. The function / (a) is known as the singularity 
spectrum of the measure ^, and we refer to /x as a multifractal measure if / (a) ^ for a continuum 
of a. In order to determine the function / (a), it was proposed to use the relationship 



/ {a) - T* (a) . 



(6) 



This relationship may not hold for a given measure. When the equality (6) is established for a measure 
/i, we say that the multifractal formalism holds for this measure. 

We recapture for geometric stationary mother process the results of Mannersalo, Norros and Riedi 
(2002), in which the conditions for £2-convergence and scaling of moments are established: 
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Theorem 2. (Mannersalo, Norros and Riedi 2002) Suppose that the conditions A' -A'" hold. 

If, for some positive numbers 5 and 7, 

exp{-5|r|} < Corr(A(0),A(.)) = ^ \Cr\-\ (7) 

then A„(t),t G [0, 1] converges in L2 if and only if 

h>l + al = M{2). 

If An{t) converges in Li, then the limit process A{t)^t G [0, 1], satisfies the recursion 

A{t) = lj\{s)dA{bs), (8) 

where the processes A (t) and A (t) are independent, and the processes A (t) and A (t) have identical 
finite- dimensional distributions. 

If Ait) is non- degenerate, the recursion (8) holds, A{1) e Cq for some q> 0, and 

00 

J2c{q,b-^)<^, (9) 

where 

c{q,t)=E sup |A«(0)-A'(s)|, 
se[o,t] 

then there exist constants C and C_ such that 

^^g-log,EA''(t) ^ E^9(i) ^ ^^9-log,EA'(t)^ (IQ) 

which will be written as 

//, on the other hand, A{1) G Cq, q > I, then the Renyi function is given by 

T{q) = q-l- log, EA" {t)^q-l- log, M{q). (11) 
If A{t) is non-degenerate, A{1) e C2, and K{t) is positively correlated, then 

VarA(t) ^ Var ( A{s)ds. (12) 
Jo 

Remark 3. The result (10) means that the process A{t), t £ [0, 1] with stationary increments behaves 
as 

\ogE [A{t + S)- A(0]' - K{q) log S + Cq (13) 
for a wide range of resolutions S with a nonlinear function 

K{q) = q - log, EA« {t)^q- log, M{q), 

where Cq is a constant. In this sense, the stochastic process A (t) is said to be multifractal. The 
function K{q), which contains the scaling parameter b and all the parameters of the marginal distri- 
bution of the stationary process process X (t), can be estimated by running the regression for a range 
of values of q. For the examples in Sections 8-14, the explicit form of K{q) is obtained. Hence these 
parameters can be estimated by minimizing the mean square error between the K(q) curve estimated 
from data and its analytical form for a range of values of q, sec Ludcna 2008, 2009. 
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3. Cq convergence: general bound 

This section contains a generalization of the basic results on multifractal products of stochastic pro- 
cesses developed in Kahane (1985, 1987) and Mannersalo, Norros and Riedi (2002). 

Consider the cumulative process A„(i) defined in (3) For fixed t, the sequence {A„(t), J^„}5^q is a 
martingale. It is well known that for q> 1^ Cq convergence is equivalent to the finiteness of 

supEA^(t) < oo. 

n 

To illustrate the approach we consider separately a simpler case 9 = 2, which was studied in Mannersalo 
et al. (2002). 
We have, 

isids2 



EAlit) = E f f A„(si)A„(s2)ds 
Jo Jo 

[|EA«(si)A(')(s2)dsids2. 



The process A^*) is stationary. Therefore 

rt ft ri 



EAlit) = 2 J]EA«(0)A«(s2-si)dsids2 

Jo Jsi^^o 

= 2 / / 'f\pib\s2-s,)dsids2 
Jo Jo ,^^0 

< 2t Y[pib'u)du, 

Jo ,_n 



where 

piu) = EA(0)A(u). (14) 
Hence, to show £2 convergence it is sufficient to show that 



pt ri 

/ Y^pih'u)d 

Jo „-_n 



sup / /9(6 M)du < CX). 

Lemma 4. Assume that piu) as defined in (14) is monotone decreasing in u, 

b > EA(0)2 (15) 

and 

00 

Y^ipiV) - 1) < cx). (16) 

i=0 

Then Anit) converges in C2 for every fixed t G [0, 1]. 

Proof. Without loss of generality let i = 1. Let n = niu) = [— logf, w] be the integer part of 
— log't, u. Then, using monotonicity of p we obtain 

n n 

n4p(&^"))<p(o)"("^ n pi^'^^y 

i=0 i=n(u) 



Using monotonicity of p again. 



n Pib'^) < I[pib''"'^^^u)) < Y[pib^u) = n. 

i=n{u) i=0 i=0 
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Constant 11 is finite due to the condition (16). For sufficiently small S £ (0, 1), by the condition (15), 
b^-^ > p(0) = EA(0)2 . Therefore, 

sup / Y[p{b'u)du<n p(0)"(")du 

n Jo Jo 

< n / < n / -^du < oo. 

~ Jo ~ Jo u'-' 

The proof of Lemma 4 is complete. 

□ 

Now we are going to consider q > 2. Here, we assume additionally that An{t) is a cadlag process. 
The starting point is the equality 



EAUt) = E / / . . . / A„(si)A„(s2) . . . K isg)d. 
Jo Jo Jo 



SidS2 ■ ■ ■ dSq 

= g! / EA„(si)A„(s2) • ■ • A„(sp)dsids2 • ■ -rfsq- (17) 

JO<Sl<...<Sg<t 

First we make change of variables 

Uo = Sl,Ul = S2 - Si, . . .Mq-l = Sq - S,_l, 

which transform equality (17) into 

PUo-i l-Mg-l <t 

E^n W = ?! / EA„(uo)A„(wo + ui) . . . A„(uo H h Uq_i)dwo . . . 

J 0<UO,...,Uq-l 



< ql EAn{uo)An{uo + ui) . . . A„{uo)duodui . . .dUq-i 
= (?! / EA„(0)A„(ui) . . . A„(iti H \- Uq^i)dui . . .duq-i, 

J 0<tti ,...,ti,_i <t 

where we used stationary of the process A{t) to obtain the latter inequality. Now let 

p(iii,...,w,_i) =EA(0)A(iii)...A(wi + ---+u,_i) (18) 

and thus it is sufficient to prove that 

« n 

sup / /9(6'mi, . . . , 6'itg_i)(iiti . . . ditq-i < oo. (19) 

We require that the function p{ui, . . . ,Uq^i) satisfies certain mixing conditions. Namely, for any 
m < q — I, 1 < ii < ii < ...,jm<'Z — 1 and Uj — 0, j ^ ii . . . , Um and Uj — A, j — ii, . . . , im, we 
assume that the following mixing condition holds: 

lim p(mi,...,u„_i)=EA(0)*iEA(0)'=-*i-...-EA(0)9-''". (20) 

A— i-oo 

We are ready now to state the main result of this section. 

Theorem 5. Suppose that conditions A'-A" hold. Assume that p{ui, . . . ,Uq^i) defined in (18) is 
monotone decreasing in all variables. Let 

b"-^ > EA(0)« (21) 

for some q > 0, and 

oo 

5](p(6",...,n-l)<oo- (22) 

n=l 

Finally assume that the mixing condition (20) holds. Then, 

EA{ty < oo, (23) 

and An{t) converges to A(t) in Lq. 
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Proof of Theorem 5. It is sufficient to prove tfiat equation (19) liolds. To simplify notation we put 
t = 1. First represent tfie integral in (19) as the sum of the integrals over different regions 

» n 
/ JJp(6'ui, . . . jfo'uq-l)^^! . . 

jO<Ul,...,«g-l<l J_Q 

^ n 

= ^ / ]^p(6'ui, . . . ,&'u,_i)(iMi . . .(i-Uq-i, (24) 

where the sum is taken over all possible permutations of numbers {1,2, . . . ,q — 1). Next we are going 
to bound the integrals on these separate regions. Put 

= Uji,M(2) = Wj^, . . ., = 

Fix a large number A > 1 which we define later and define an auxiliary function n{u) = — [logt, u/A]. 
Note that this function is no n- negative for u < 1. Now let 

/i = n(u(i)), h = n-(w(2)), • ■ • , lq~i = ?^(^i(<^-l))• 
These numbers are decreasing 

^1 > ^2 > •■• > ^9-1. (25) 

Then we can split the product as 

i[pib'ui,...,b'u,.i)= n n ■•• n n^(^'"i'---'^'"9-i)- (^6) 

1 = 1 = l = lq-l l = h l = li 

Further, using monotonicity of the function p wc can estimate for I < Iq-i, 

p{h^ui, b^Ug_i) < p{0, . . . , 0) = EA(O)'?. 

For I G [lq^i,lq-2), we have 

p{b'ui,...,b'uq_i) < p{0,...,0,A,0...,0), 

where iq_ith argument of the function p is equal to A and all other arguments arc equal to 0. Indeed 
this holds due to the fact that for / > Iq^i 

b'^U^q_^ > > "(g-l) = A 

and the monotonicity of the function p. Here recall that U(<j-i) corresponds to Ui^ -^. Fix a small 
number 5 which wc define later. Now wc can note that mixing condition (20) implies that 

lim p(0, . . . , 0, A, 0, . . . , 0) = EA(0)''-iEA(0)«-*'-i 

A^oo 

Hence we can pick A = A{5) sufficiently large to ensure that 

p{{),...,Q,A,{),...,Q) < (l + (5)EA(0)*'-iEA(0)''-'''-^ 

Function g{x) = InEA(O)^ is convex. Hence we can apply Karamata majorization inequality [32] to 
obtain that 

5(i«-i) + 9{q - ig-i) < diq - 1) + 3(1)- 

Therefore, 

EA(0)*'-iEA(0)''-'''-^ < EA(0)'?-iEA(0) = EA(0)«^^ 

and 

p(0,...,0,A,0,...,0) < (l + (5)EA(0)«-i. 
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Similarly, for I S [lq^2,lq-3], we have 

p{b'ui,...,b'uq^i) < p(0,...,0,A,...,0,A,0...,0), 

where the arguments of the function p are equal to except arguments iq-i and iq-2 which are equal to 
A. Applying the mixing condition and increasing A if necessary we can ensure that for I G [lq-2, Iq-s], 

p{b'ui, b^uq-i) < (1 + 5)EA(0)"EA(0)^-'^EA(0)'^-^ 

where a = min(ig_2, ig^i), b — max(iq_2, ig-i)- Wc apply now Karamata's majorisation inequality 
twice. First application of the inequality gives 

EA(0)''EA(0)^-" < EA(0)^-^ 

Second application of Karamata's inequality gives 

EA(0)^-iEA(0)«-'' < EA(0)9-^ 

Hence, for I G [lq~2,lq-3) and sufficiently large A, 

p{b^ui, . . . , b^u"-^) < (1 + (5)EA(0)''EA(0)''-"EA(0)9-'' < (1 + 5)EA(0)«-^ 

In exactly the same manner, using the mixing conditions and Karamata's majorisation inequality one 
can obtain for / £ [Ij, Ij-i) and j = q — l,q — 2, . . . ,2 

p{b^ui, . . . , b^ui-^) < (1 + ^)EA(Oy . 

Hence, 

h-l i<,_2-l 

np(&'ui, n n •••n^(^'"i'---'^'"«-i) (2?) 

1 = 1 = l = lq^l l = l2 

q 'i-l — 1 q 

<{l + 5)''X{ n EA(0r = (l + (5)'i J];(EA(O)0''"'"'\ 

i=2 l=li i=2 

where Iq — Q. Rearranging the terms we can represent this product in a slightly different form 
Now one can note that since U are decreasing, see (25), 

Iq-l + ---+k< ^^—^{ll + ■■■+ Iq-l), 

for any j = — 1. Indeed, the latter inequality is equivalent to 

{i - l){lq^i + ---+h)<{q- i){h-i + ■■■ + h), 

which follows from 

Iq-l + ' ' ' + ^ ; ^ ; ^ + ' ' ' + ^1 

: S H S S : • 

q ~ t I — 1 

In addition, by the Karamata's majorization inequality, 

EA(0)'EA(0)'-2 



(EA(0)*-i)2 

Therefore, 



> 1. 



2\ '<j-i+---+'' /'p.\(n\iv.MnY-2\ 



I^EA(0)*EA(0)*-2y<'-i+-+'- ^ ^EA(0)'EA(0)'' 



V (EA(0)'-i)2 ; - V (EA(0)'-i)2 
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Hence we can continue (28) as follows 



n(EA(or)'--''-^<n 



i=2 



\ (EA(0)'-i)2 

. , q-i + l-2(q-i) + q-i-l ij + .-. + l 

(EA(0)) — — Jl (EA(O)'-i) '-^ = (EA(0)) — — . (29) 



i=2 

Plugging the latter estimate in (27) we arrive at 

h 



J{p{b'ui,...Mu^-i) < {1 + 5)'' (EA(0))-'^' 



We can now make use of the condition (21) and by taking 5 sufficiently small we can ensure that 

h 

W p{b'ui, . . . , b'Ug.l) < bii-^)ih + - + lq-l) ^ (^^^2 . . . W,_l)-1 + M«(l-") (30) 

for some small £ > 0. We are left to estimate the product 11"=;! uniformly in n. For that we are going 
to use finiteness of the series in (22). First note that for I > h, 



Uu,j > Ah' 



i-h 



Then, by monotonicity of the function p, uniformly in n, for some C > 

n n 

l=ll l^h 

OO 

<]lp{b\...,...,b')<C, (31) 

1=0 

according to the finiteness of the series. Together (30) and (31) give us 



. n 

/ Y[p{b'ui,...,b'ug. 

J0<«1,...,M,-1<1 J_Q 



i)dui . . .du 



9-1 



» n 

= ^ / J^p(6'mi, . . . , 6'ug_i)cfui . . . duq-i 

< C ^ / {uiU2 ■ ■ .Uq^l)~^^''dui . . .dUq-l, 

il,. "'0<«ii<U,2<U,^_j<l 

= C {uiU2 ■ ■ ■Uq^l)~^^^dui . . .dUq^l (32) 

which immediately gives a finite bound for EA„(1)* uniform in n. 

□ 

Remark 6. It is not difficult to show that (21) is sharp. Indeed suppose that 

6«"i < EA(0)« 
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and that p(ui, . . . , is continuous at (0, . . . , 0). Then, for e > 0, 

/•Uq-\ hMq-l <t 

EAlit) ^ql EA„(0)A„(iti) . . . A„(mi + • • • + u^-^duo . . . dug^i 

'JO<UQ,...,Uq^l 

fua-l <t " 



/•Moi h«5-l 

= g! / . . . ,6'ug_i)(iwo ■ • -c^Mg-i 

0<Uo,...,«5-l ;_Q 

> ql ]^p(6'mi, . . . ,6'-Uq_i)(i-uo • ■ -dw, 

J0<Mo<l/2,0<Mi...,ti,_i<e/&" 



> ^ / J]^ • . ■ , e)dtti . . . dug^i 

9! 9-1 f pie, ■■■,£] 



2 V 

Since p(e, . . . , e) can be made arbitrarily close to p{0, . . . , 0) — EA(0)^, then, for sufficiently small 
e > 0, 

p{e,...,e)>b''-' 

and 

as n — )• cxD. 

4. Scaling of moments 

The aim of this Section is to establish the scaling property (10). For q > 1 let 

, /EA(0)«-iA(sit) \ , , 

Note that pq{s) < 0. For g e (0, 1) let 

, , /EA(0)9-iA(su) \ , ^, 

For (J < 1 it is easy to see that Pq{s) > 0. 

Theorem 7. Assume that A{t) G £g and Pq{s) defined in (33) and (34) is such that 

00 

J2 Mb-n\ < 00. (35) 

Then, 

EA« (t) - t«->°S'' t e [0, 1] . (36) 

Proof of Theorem 7 

Our strategy in proving of (36) is to use martingale properties of the sequence An{t). We concentrate 
mainly on 5 > 1, as the case g < 1 is symmetric. For the upper bound we obtain uniform in n 
estimates from above for EAn{t)'^. Then, since A„(t) converges to A{t) in Cq, the same estimates hold 
for EA(t)'. For the lower bound, wc use the fact that as An{t) S £g for g > 1 the martingale A„(t) 
is closable. Hence it can be represented as An{t) = E{A{t)\Ai{t), . . . , A„(t)). Therefore, for q > 1, by 
the conditional Jensen inequality, 

EAnit)" - E(E(^(i)|Ai(t), . . . , A^it)))" 

< E{E{A{tf\Ai{t), . . . , Anit))) = EA{ty. 
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Thus, we are going to obtain an estimate from below for F,An{t)'' for a suitable choice of n. Clearly, 
by the latter inequality, this estimate will hold for EA(i)'? as well. 
We start with a change of variable 

An{t) = f A„(s)ds 

JO 

= t An{ut)du = tA,i{t). 
Jo 

Clearly An{t) is a martingale for any fixed t. 

We are going to treat the cases q > 1 and g < 1 separately. This is due to the fact that for 
q > 1, the sequences An{t)'^ and An(t)'^ are submartingales while for q G (0,1) the sequences are 
supermartingales with respect to the filtration J^n = o-{A^^\ . . . , A'"'). 

We start with an upper bound for 5 > 1. Let nt = — [log^i] be the biggest integer such that 
Tit < — log;, t. We use the Holder inequality in the form, 

"l/»l)' = (/l/IM"«M"f<(/"l/™) (/h)'". 

where 1/q + 1/p = 1. It follows from the latter inequality, 

in 



k=0 / 

(^n^*"("^)j l[ ^^''\^t)duj i^j l[ A^''\ut)duj . 
Applying expectation to both sides we obtain, using independence of A^^^ of each other, 

„1 nt — 1 n / / pi n ^ 



/ \[ E(A('^-))«(wt) \[ E A'^''\ut) / \[ A^^\vt)dv 

"'O k=0 k=nt \ V" k=nt / 

By the stationarity of the process A(t) we have 

nt — 1 



du 



J]^ E(AW)9(wt) = (EA(O)'?)"' < (EA(0)«)"'°*^'* ^ ^-iog,EA(o)' 

k=0 



Therefore, 



,1 n \ 



/in \ 

^^-iog,EA(o)^E j^y [] AW(^.i)rf«j 

/ »i n-nt \ 9 



fc=0 
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Now note that 



< & sup EA(s)«. 

se[o,i] 

This bound is uniform in n and therefore, 

se[04] 

Now we turn to the lower bound for q > 1. 
Since A„(t) is a submartingale, 

El(t)«>El„,(t)«, 

where rif ~ [— log;, . 

We are going to obtain a recursive estimate for EA„(i). First, 

= E ^ A„(ut)(A("+i' (6"+^ut) - A("+i) (0))dw + I„(0 A("+i) (0)^ 
Now we can use an elementary estimate of the form: if a + 6 > and b > then 



(37) 



for q > 1. This estimate is easy to prove by analyzing the function (l + t)'' — 1 — qt for t > — 1. Applying 
(37) we obtain 



9-1 



EA„+i(t)9 > gE 

+ E(l„(t)A("+i)(0)j^ =£i+S2. 
The second expectation is straightforward. 



A„(t)A("+i)(0)) / A„(Mi)(A("+i)(&"+iut)- A("+i)(0))du 



-E2 — E 



A„(ut)A("+i)(0)dM I = EA(0)«EA„(t)«, 



(38) 



(39) 



where we use independence of A„ and A("+^^. For the first expectation, rearranging the terms, we 
have 



Ei^qE I ^„(i)*-iA„(wt)(A("+i)(0))«-i(A("+i)(6"+iui)- A("+i)(0))du 

El„(t)«-iA„(wt)E(A("+^)(0))«-nA("+iH&"+iut) - A("+i)(0))du 
By the definition of pq, see (33), for all u E [0, 1], 

E(A("+i)(0))9-i(A("+i)(6"+iut) - A("+i)(0)) > EA(0)Vg(^^"+^0- 

Therefore, 

Ei>qf EAn{t)''-^An{ut)duEA{oypg{b"+H) 



> qE 



A„ (<)'-! / A„(tit)d^ 



EA(0)V9(&"+'t) 



gEA„(0'EA(0)V9(&"+H) 
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Therefore 

Si >9EA„(0«EA(0)V,(&"""')- 
The latter inequality together with (38) and (39) gives us 

El„+i(t)« > El„(t)''EA(0)« (1 + gp,(6"-"')) (40) 

Now we can iterate it. First fix N* such that \qpq{b^^^)\ < I for n > N*. Then, iterating (40), we 
obtain 

nt-N* 

E<_^. > (EA(0)')"-^* n (l + gp,(6"-"')) 

oo 

> (EA(0)')"-^* J] (f + gp,(6-")) 

ft is sufficient to note that the latter product is strictly positive due to (35). As A„(t)5 is a submartin- 
gale, we have 'EA{t)'^ > EA*^_^. and the required lower bound for g > f follows. 

The proof for q € (0, 1) is symmetric. For these values of q and a fixed t, the process An{t)'^ is a 
supermartingale with respect to the natural filtration J^n ~ o'(A'-^-', . . . , A*^"'). The bound from below 
is proved using the reverse Holder inequality for q G (0, 1) and p such that 1/p + 1/q = 1. 



1/9 / ,\ \ 1/P 



Note that p is negative. We are going to use this inequality in the form, 

' = (J^ > i/n.9i) ( t \9\ 







It follows from the latter inequality, 

' nl " 



]J Afc(ut)di 



fc=0 



> 



V"'" \fe=0 / fe=nt+l / \"^° k=nt + \ 



The rest of the proof goes exactly as the proof of the upper bound for g > 1. 

To prove the upper bound, we proceed similarly to the proof of the lower bound for q> \. First we 
establish a recursive estimate. The elementary inequality (37) still holds (in the opposite direction), 
forge (0,1), 

(a + hf < qab'^-^ + b" 

for a + b > 0,b > 0. Repeating step by step the arguments for g > 1 we obtain an upper bound 

eI„+i(<)« < El„(t)«EA(0)' (1 + 9/5,(5"-"*)) . 
Applying this bound recursively 

^Al_^, < (EA(0)')"*-^* n (1 + '?P9(^"""')) 

n=0 

< (EA(0)')"-^* f[ {l + qp,{b-)) 

It is sufficient to note that the latter product converge due to (35). As An(t)'' is a supermartingale, 
we have F,A{t)'^ < EA'^^ and the required upper bound for q < I follows. 

□ 
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5. Log-normal scenario with possible long-range dependence 

The log-normal hypothesis of Kolniogorov [34] features prominently in turbulent cascades. In this 
section, we provide some related model, namely the log-normal scenario, for multifractal products of 
stochastic processes. In fact, this log-normal scenario has its origin in Kahane [30, 31], see also [1]. 
In this section we present a most general result on log-normal scenario for a models with possible 
long-range dependence. 

We combine Theorems 5 and 7 for this special case in order to have a precise scaling law for the 
moments. 

B'. Consider a mother process of the form 

A(i)=exp(x(t)-ia2,\, (41) 



where X(t), t G [0, 1] is a zero-mean Gaussian, measurable, separable stochastic process with covari- 
ance function 

Rx{T)^aj,CoTr{X(t),X{t + T))) 

Under condition B', we obtain the following specifications of the moment generating functions (1) and 
(2): 



M(C) = Eexp<{c(X(t)--a| ) ^ =e^--(C--C)^ 



1 

— ( 
2 

Af(Ci,C2;ii-t2) = Ecxpjci (^X(ti)-lai^ +C2 (^^(t2)-^^l: 

= expjidi [C?-Ci+C|-C2] +CiC2i?Jf(ii-t2)|, Ci,C2eR\ 

where cr^ e (0, oo). It turns out that, in this case, 

M(l) = 1; Af(2)=e'"i; al ^ e"^^ ~ 1; 
Cov(A(ii),A(t2)) = M(l,l;ii-t2)-l = e«^(*i-*^)-l 



and 



log,EA(i)« = i<-^, q>0. 
2 log 



Note that 

gflx(ti-t2) _ I > ^^(^^ _ ^2). 

Using Theorem 2. we obtain 
Theorem 8. Suppose that condition B' holds with the correlation function 

CoTT{X{t),X{t + T)) <Ct-°', a>0, (42) 

for sufficiently large r, and for some a > 0, 

l-CorT{X{t),X{t + T)) KCItI" , (43) 

for sufficiently small r. Assume that 

6>exp{gVi/2}, (44) 

where q* > is a fixed integer. Then the stochastic processes 

rt " 



r J. 

•^0 j=o 
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converge in Cq, < q < q* to the stochastic process A{t), t G [0, 1], as n — > oo, such that 

^A{tY - t-'^'?'+(''+l)9^ q e [0, q*l (45) 
and the Renyi function is given by 

T{q) = -aq' + ia + l)q-l,qeiO,q*), 

where 



2 log 6 

Moreover, if 

Corr(X(t), X{t + r)) = > 0, 

where L is a slowly varying at infinity function, bounded on every bounded interval, then 

o ^ 1 L(t\u — v\)dudw „ ^ ,,„, 

VarA{t)^t^~''ol / \\ ,0<«<1, 46 

Jo Jo L{t)\u-w\ 



VarAit) ^ 2tal / (1 - I)^^dr,a > 1. 
Jo £ V\ 



(47) 



Remark 9. We interpret the inequality (46) as a form of long-range dependence of the limiting process 
in the following sense: one can replace the interval [0,1] into more general interval [0,t], and for a 
large t we have the following: 

.. . /■* f*' L{\u — v\)dudw ,. o 9 L(t\u — vDdudw 

VarA(t) > lim / / = lim t^-'^aj^ / / v I 1/ 



i-s-oojQ Jq |u — w| t->oo Jq Jq L{t) \u — w\ ' 
and 

L(t\u~v\)dudw dudw 2 

lim / / ^TTTTI ,a = / / 1 [a =7^ T77^ -,0<a<l, 



*^°°7o 7o i(t)|u-wr Jo Jo \u-w\°' (l-Q)(2-a) 
while the inequality (47) has the form of short-range dependence of the limit process: 

pt ^ L(t) f°° L(t) 

lim / (1 ) I ,a dr ~ / , dr < oo,a > 1. 

t^oo_/o t \t\ Jo \t\ 

For a ~ 1, one can show that for some other slowly varying at infinity function Li[t), bounded on 
every bounded interval, 

tUt) 



lim 

t^oo 



/ (1--)^ 

io^ \r\ 



■dr 



lL,{t) = 1. 



Remark 10. Note that the correlation function CoTT{X{t), X{t + t)) = (1 + |r| ) a > 0, satisfies 
all assumptions of the Theorem 2, among the others. 

Proof. To prove Cq convergence wc use Theorem 5. The moment generating function of the multidi- 
mensional normal distribution is given by the following expression 

M(Ci, C2, . . . , O = Ee';^^(^^)+-+«''^(^') = cxp I icRC*} , 

where 

C = (Ci,-.-,C,) 

and R is the covariance matrix 

R={i?x(|s«-Sj|)}»j=i..9 



D. Denisov and N. Leonenko/ Limit Theorems for Multifractal Products 



16 



One can immediately see that 

E(A(.si)A(.S2)...A(sJ) = Ee^('i)-3'^i ...e^("')-2''i 



M(l,l,...,l)e"^'^x 

e^ELi 1:1=1 «x{|3,-^,l)g-f<Ti 



We can now substitute this into ( 18) and obtain 



p{ui,U2, ■ ■ ■ ,Uq_l) = exp ^ RxiUi-\ h Uj) 

[l<l<3<q-l 

Since the function Rx{u) is monotone decreasing in u, function . . . , Uq^i) is monotone decreasing 
in all arguments. Next we need to check the mixing condition (20). Let 1 < ii < i2 < ■ ■ • < *m and 
Ui = A if i Cz {ii, . . . , im} and otherwise. Then, as ^ — >■ 00, and iq = 0, i,„+i ~ q 



lim p{ui, . . . , Uq-i) = exp 

l<k<m+l ik-l<i<j<ik 



EA(0)*iEA(0)^^-'i • ...EA(0)''-^"\ (48) 



where we used that EA(0)' e ' 2 ''^x. Finally, we should check the convergence of the series (22). 



exp{i?x (<?&"} < p(6", ...,&")< exp ^ 



We have, 



As n -> 00, Rx{V^) 0. Hence 

p{b\ &") - 1 < (1 + o(l))fcili?^(&") 

and 

p(6",...,6")-l> (l + o(l))i?x(g6"). 

As both sums 

00 00 

Rxiqb") < 00, J2 Rxib") < 00, 

n—l n—1 

the convergence of the series (22) follows. Condition (21) becomes 

-1 J^^/n^Q „™ /'7('7- 1) 2 



>EA(0)«=exp|i^Y^ai 

which is equivalent to (44). 

Next we are going to prove scaling (45). For that we apply the results of Section 4. We need to 
show that (35) holds for pg, where q £ {0,q*) and pq is defined in (33) and (34). For q > 1 we have, 
for sufficiently small s, 

\pqis)\ = - inf fEA(0)^-U(,s.) \_ r.^^a,-i)P.isu)^i-,, _ A 
"^'^ u<i \ EA(0)9 J u<i\ J 

< sup fl - e(i-9)'^x(^«)"') < 1 _ e{i-g)'^li'^r < _ l)<Tis^ 

u<l ^ ' 

Thus using condition (43) one can immediately see that the series (35) converges. For 5 < 1, the same 
arguments give bound 

Pq[s) < il-q)aj,s''. 

Using condition (43) one can immediately see that the series (35) converges. Therefore, by the results 
of Section 4 scaling (45) holds. □ 
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6. Geometric Ornstein-Uhlenbeck processes 

This section reviews a number of known results on Levy processes (see Skorokhod 1991, Bertoin 1996, 
Kyprianou 2006) and OU type processes (see Barndorff- Nielsen 2001, Barndorff-Nielsen and Shephard 
2001) The geometric OU type processes have been studied also by Matsui and Shieh 2009. 
As standard notation we will write 

k(z) = C{z;X} = logEexp{izX}, zeM 

for the cumulant function of a random variable X, and 

A'{C;X} ==logEexp{CX}, C e D C C 

for the Levy exponent or Laplace transform or cumulant generating function of the random variable 
X. Its domain D includes the imaginary axis and frequently larger areas. 

A random variable X is infinitely divisible if its cumulant function has the Levy-Khintchine form 

C {z;X} = iaz- -z^ + I (e*^" - 1 - izul[_i^i] (u)) (dw) , (49) 

where a £ M, d > and v is the Levy measure, that is, a non-negative measure on M such that 

i/({0})=0, j mm{l,u^)i'{du) <oo. (50) 

The triplet (a, d, v) uniquely determines the random variable X. For a Gaussian random variable 
X ~ iV (a, d) , the Levy triplet takes the form (a, d, 0) . 

A random variable X is self-decomposable if, for all c £ (0, 1) , the characteristic function / (z) of X 
can be factorized as f (z) ^ f (cz) fc (z) for some characteristic function fc (z) , z G R. A homogeneous 
Levy process Z = {Z (t) , t > 0} is a continuous (in probability), cadlag process with independent 
and stationary increments and Z (0) = (recalling that a cadlag process has right-continuous sample 
paths with existing left limits.) For such processes we have C {z; Z (t)} ~ tC {z; Z (1)} and Z (1) has 
the Levy-Khintchine representation (49). 

If X is self-decomposable, then there exists a stationary stochastic process {X (t) , t > 0}, such 

that X{t) =X and 

X{t)^e-^'X{0)+ f e-^^'-'^dZ{Xs) (51) 

J{0,t] 

for all A > (see Barndorff-Nielsen 1998). Conversely, if {X (t) ,t> 0} is a stationary process and 
{Z{t),t>0} is a Levy process, independent of X {0) , such that X (t) and Z (t) satisfy the Ito 
stochastic differential equation 

dX (t) = -XX (t) dt + dZ (Xt) (52) 

for all A > 0, then X (t) is self-decomposable. A stationary process X (t) of this kind is said to be an OU 
type process. The process Z (t) is termed the background driving Levy process (BDLP) corresponding 
to the process X (t) . In fact (51) is the unique (up to indistinguishability) strong solution to Eq. (52) 
(Sato 1999, Section 17). The meaning of the stochastic integral in (51) was detailed in Applebaum 
(2009, p. 214). 

Let X (t) be a square integrable OU process. Then X (t) has the correlation function 

Corr(A:(0), A:(t)) = rx (t) = exp {-A \t\} . (53) 
The cumulant transforms oi X = X{t) and Z (1) are related by 

C{z;X}= r C{er^z-Z{l)]ds= f C{i-Z{l)}'^ 
Jo Jo s 

and 

dC{z;X} 
C{z;Z{l)} = z — . 
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Suppose that the Levy measure v oi X has a density function p (u) , m S M, which is differentiable. 
Then the Levy measure of Z (1) has a density function q (u) , u € K, and pand q are related by 

q (u) = —p (u) — up' (u) (54) 

(see Barndorff- Nielsen 1998). 

The logarithm of the characteristic function of a random vector {X{ti), ...,X{tm)) is of the form 

\ogEcxp{i{ziX{ti) + ... + z,nX{t^)}^ / K(^Zj-e-^(*^-")l[o,oo)(ij (55) 

where 

k(z) = logEcxp{izZ(l)} ^C{z;Z{l)}, 

and the function (55) has the form (49) with Levy triplet (a, d, v) of Z{1). 

The logarithms of the moment generation functions (if they exist) take the forms 

logEcxp{CX(t)} = Ca + -C'+ / (e«"-l-Cul[-i,i] {u))^ [du) , 



2 

where the triplet (a, d, z^) is the Levy triplet of ^(0), or in terms of the Levy triplet (a, d, v) of Z{1) 
logEexp{CX(t)} = a^(Ce-^(*-^)l[o,„,)(t - s))ds + ^C' ^(Ce-^'*-^)l[o,oo)(< - s))^^^ 
+ / /[cxp|Ke-^(*-^)l[o,„,)(t-s)|-l-^.(Ce"^(*-^)l[o,^)(i-s))l[_i,i](M)]i>(d7.)ds, (56) 



and 



logEcxp{CiX(ti) + C2^(t2)} 

= « / (EOe-^'*^"^'l[o,oo)(t,-s))ds + |c' / (^Oe-'^*^-^)l[o.oo)(i,-s))') 



[cxpL;^ge"^'*^"'^l[o,oo)(tj -s) i - 1 



2 



^Oe-^(*^-^'l[o,oo)(i,-s)j l[_i.i](7.)]j>(d7.)ds. (57) 

Let us consider a geometric OU-type process as the mother process: 

K{t) = e-^f'^-'^^.cx = logEe^(°\Af(C) = Ee^^^^*)-'^-^), Afo(C) = Ee'^'^^*) 
where X{t),t e M+, is the OU-typc stationary process (51). Note that 

Mojq) ^ M{q) 
Afo(l)9 M{l)i- 
Then the correlation function of the mother process is of the form. 

Corr(A(i), A(i + r)) = ^^[^'^^'^^^ ~ \ (58) 

where now 

M(Ci, C2; t) = Ecxp{Ci(X(ti) - cx) + UX{t2) - cx)} 

= exp {-(Ci + C2)cx} Eexp{CiX(ti) + C2^(i2)}, (59) 

and Ecxp{CiX(ti) + C2^(t2)} is defined by (57). 

To prove that a geometric OU process satisfies the covariance decay condition (35) in Theorem 7, 
the expression given by (57) is not ready to yield the decay as i2 — — > 00. 

The following result plays a key role in multifractal analysis of geometric OU processes. 
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Theorem 11. Let X(t),t G M+ be an OU-type stationary process (51) such that the Levy measure v 
in (4-9) of the random variable X(0) satisfies the condition: for a positive integer q* € N, 



Then, for any fixed b such that 



xei ""vi^dx) < oo. (60) 

x\>l 



Mo(l)'3 



stochastic processes 

An{t) = / n^''^ e [0,1] 

converge in Cq to the stochastic process A[t) € £g, as n ^ oo, for every fixed t G [0, 1]. The limiting 
process A{t),t g [0, 1] satisfies 

EA«(t) - t«-'°S''^'^'(*), ge[0,g*]. 

The Renyi function is given by 

T{q) = g - 1 - logfc EA« {t) = q[l + - log, Mo{q) - 1, g G [0, q*]. (62) 

In addition, 

VaTA{t)^2tJ^ (l-^) (Af(l,l;T)-l)dr, (63) 

where the bivariate moment generating function M(^i,(^2',ti — t2) is given by (59) 

Proof of Theorem 11 We are starting with Cq convergence. To show the convergence we apply 
Theorem 5. It is sufficient to show the convergence for q = q* since the convergence for q < q* 
immediately follows from the convergence for q = q* . First we will derive a suitable explicit expression 

for p{ui, . . . , Uq-l). Put Sl = < S2 ~ Ul < S2 = Ul + U2, ■ ■ ■ , Sq = Ul + ■ ■ ■ + ThcU, 

p{ui, Uq^i) = EA(si) . . . A{sq) = E cxp {X{si) + ... + X{sq) - qcx} ■ 
Using representation (51) one can obtain 

X{sq) = e-^^'-'-'-'-^^Xisq-i) + [ e-^^'-'-'UziXs). 



'(s,„i,s,] 

Then, using independence of X{sq^i) and the integral /^^ ^ ^ j e~^^*''~*^(i2'(As) we obtain 

Eexp{X(si) + ...+X(s,)} 

= Eexp {X(si) + . . . + (1 + e-^("'-"'-i))X(.s,_i)} Ee^<=«-i-<'i ^"'''''"°'''^(^^) 

= Eexp {Xis,) + . . . + (1 + e-(----))X(.,_,)} ^^jf^,,,^ 
= Eexp {Xis,) + ... + (! + e-(— ))x(.,_,)} ■ 

Proceeding further by induction we obtain 
Eexp{X(si) + ... + X(s,)} 

- M (1) ^ e-^"'-i)Mo(l + 6-^""-= + e-^("''-i+"'-=)) . . . Mo(l + e"^"i + . . . + e-^("i+-+"5-i)) 
~ A'/o(e"^"''-OA/o(e"^"'-" + e--^("''-i+"<'-2) . . . Afo(e"^"i + . . . + e-^("i+---+"'-i)) 
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Hence 



p(ui, . . . ,Ug_i) = 



Afo(l)Mo(e-^"i+-+^"''"'+-+"''-'') 

This representation allows us to show monotonicity of p{ui, . . . ,Uq-i). For that we use the following 
inequality 

W + < Mo{l + t) 
Mo{s) - Mo{t) ^^^> 

for s < t. This inequality follows from the fact that lnA/o(i) is a convex function and the Karamata 
majorisation inequality. Hence 

A/o(l + .s) 
Mo(l)A/o(s) 

is monotone increasing in s. Since e~^'^ is monotone decreasing in u the representation (64) implies 
that p(ui, Uq-i) is monotone decreasing in all variables. 
Condition (21) becomes 

7-1 . T.A/n^ _ Moiq) 



b"-' > EA(0) = 



Mo(l)«' 

which is equivalent to (61). 

To show the finitencss of the series (22) we are going to use the following statement. 

Lemma 12. For s G [0, 1], the following estimate holds 



Mo(l)M(s) - VAfo(l)eE^(i) 
Proof of Lemma 12 Function InAio(i) is convex. Therefore, 

lnA/o(l + s) = lnA/o((l - s) + 2s) < (1 - s)lnA/o(l) + slnA/o(2) 
In addition, by the Jensen inequality, 

Mo(s)=Ee^^« >e^E^(i). 

Together these inequalities imply, 

Mo(l + s) ^ A/o(l)i-^A/o(2)^ _ ( A/o(2) 



Aro(l)Mo(s) - Afo(l)e^EX(i) \M^{l)e^x(i 



Now, using (64) and monotone decrease of A'/o(l + s)/Afo(s) 



□ 



^- U^o(l)Afo(e-^^")y ^ ^ 

< C""" < l + o(l)lnCge" 



-A6" 



where the former inequality follows from Lemma 12 with C = A/o(2)/(A/o(l)e^^^^^). Then, conver- 
gence of the series (64) follows from the finiteness of the series X^^^i ^^^^ ■ 

Finally we need to check the mixing condition (20). Let 1 < ii < 12 < ■ • • < *m and Ui = A ii 
i G {ii, . . . , im} and otherwise. In this context it is convenient to use (64) in the form 



p{ui,...,Uq-i) = Y[ 



/j^iAfo(l)A/o(EL^"'^'-"') 
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Then, as oo, and io = 0, i,n+i = q 



^o(l) 



This proves (20). Therefore Theorem 5 gives £q convergence of An{t). 

To prove the scahng property we are going to use the results of Theorem 7. First using representation 
(51), we have for any q, 

EA(t) A(0)«"i = E exp{(q - 1)A:(0) + X{t)- qcx} 

= Eexp{(q - 1 + e-^*)X(0) + / e-^^'~'Uz{\s) - qcx} 

Eexp{e-^*X(0)+/(o_,)e-^(*-«)dZ(As)} 



Ecxp{(q - 1 + e-^')X(O) - qcx} g^F^^ 

Eexp{(,-l + e--)X(0)-,c,}|pg|i 
Mo(g - 1 + e-^*) 



Then, 



Mo(l)9-iM(e-^*)' 

EA(t)9-iA(0)9-i Ma{q - 1 + e^^*)Mo{l) 



EA(0)9 Afo(g)Mo(e-^*) 

For q > 1 the latter function is monotone decreasing in t, as follows from the Karamata motorization 
inequality. Hence, 

^ I — my^j ^ ' — Mo(,)Mo(e--) • 

Function /(x) = lnM(a;) is convex. Condition (60) ensures that the derivative /'(<?) exists for q < q* . 
Then, for any x < q, 

f{x)-f{q)>{x-q)f'{q). 

In particular for x = q — 1 + e""^*, 

f{q-l + e-^n - fiq) > i-l + e-^^)/'(g). 

In addition, by the Jensen inequality, 

Mo(e-^^) = Ee^""^(°) < (Ee^(°))^"" = A/o(l)^"". 

The latter two inequalities give 

\Pgis)\ < l-e(-i+^^'°)(/'(9)-/(i)) (70) 
< (1 - e-^-)ifiq) /(I)) < Xsifiq) - /(I)). 

Then 

OO OO 

< J2 \P^(b~")\ ^ M/'(<?) - /(I)) E < ^- 

n—l n—1 

Since we have already shown that A{t) G Cq we can apply Theorem 7. 

□ 
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7. Superpositions of geometric Ornestein-Uhlenbeck processes 

The correlation structures found in applications may be more complex than the exponential decreasing 
autocorrelation of the form (53). Barndorff- Nielsen (1998) (see also Barndorff-Nielsen and Sheppard 
(2001)) proposed to consider the following class of autocovariance functions: 

m 

i?™(i) -^fT|cxp{-A, (71) 

which is flexible and can be fitted to many autocovariance functions arising in applications. 

In order to obtain models with dependence structure (71) and given marginal density with finite 
variance, we consider stochastic processes defined by 

dXj (t) = -XjXj (t) dt + dZj (Xjt) , j = 1, 2, m, ... 

and their finite superposition 

^msup(t) = Xi{t) + ... + Xrn(t), t > 0, (72) 

where Zj, j = 1, 2, m, ... are mutually independent Levy processes. Then the solution Xj = 
{Xj(t),t > 0} , j = 1,2, ...,m, is a stationary process. Its correlation function is of the exponential 
form (assuming finite variance of the components). 

The superposition (72) has its marginal density given by that of the random variable 

^™sup(0) = Xi(0) + ... + X,„(0), (73) 

and autocovariance function (71). One can generalize the Theorem 10 to the case of finite superposition 
process (72). 

We are interested in the case when the distribution of (73) is tractable, for instance when X„i{0) 
belongs to the same class as Xj{0),j = 1, m (see the examples in Sections 8-13 below). 
Define the mother process as geometric process 

A(t) = e-^"-p(*)"'=^,cx = logEe^"™p(°\ A/(C) = Ee'^(^'"="p(*)-'=^\ Mo(C) = Ee'^'^'"="p(*), 

where Xmsup{t),t G K+, is the finite superposition process (72). Note that 

m 

logEexp{ClX„,up(ti) + C2Xmsup(t2)} = ^ log E CXp{ClXj (^i) + C2Xj{t2)}, 

J = l 

where logEexp{CiXj (ii) + ^2-'^^ (^2)}, .? = 1, ••, m are given by (57). 
Denote 

M(Cl,C2;tl - t2) = exp{-Cx(Cl + C2)}Ecxp{ClX„,up(ti) + C2^msup(i2)} (74) 

We can formulate the following theorem which can be proved similar to Theorem 11. 

Theorem 13. Let Xmsupit),t € IR+ be a finite superposition of OU-type stationary processes (72) 
such that the Levy measure v in (49) of the random variable Xm{t) satisfies the condition that for a 
positive integer (7* £ N, 

/ xe'^'''iy{dx) < 00. (75) 

J\x\>l 

Then, for any fixed b such that 



stochastic processes 



.Mo(l)«* 
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converge in Cq to the stochastic process A{t) G Cq, as n —> oo. The limiting process A{t) satisfies 

EA9(i) qe [0,q*],te [0,1]. 

The Renyi function is given by 



T{q)^q-l-\og,EAHt)=q{l + :^)-log,Mo{q)~l, qe[0,q*]. (77) 

log 6 



In addition, 



VarA(t)^2t^ (l - ^) A/(Ci, C2; (78) 

where the bivariate moment generating function M(Ci,C2;ii — t2) given by (74) 

We are interesting to generalization of the above result to the case of infinite superposition of 
OU-type processes which has a long-range dependence property. 

Note that an infinite superposition (m — > 00) gives a complete monotone class of covariance func- 
tions 

/•OO 

i?sup(i) - / e-*"rfC/(u),i>0, 



for some finite measure U, which display long-range dependence (see Barndorff-Nielsen 1998, 2001, 
Barndorff-Nielsen and Leonenko 2005 for possible covariance structures and spectral densities and 
Barndorff-Nielsen, Eiler and Stelzer, R., 2011 for multivariate generalizations). 

We are going to consider an infinite superposition of the OU processes, which corresponds to 
m = 00, that is now 



assuming that 



In this case 



EXj (t) <oo,Y^ YarXj (t) < 00, (80) 
3=1 i=i 

00 

Rsnp{t)^J2'^h^Pi-^^\*\}^ (81) 



and if we assume that for some 5j > 



EXj{t) = S,Ci,VaTX,{t) = g] = = ^-d+^d-^)), i < i/ < 1, (82) 

z 

where the constants Ci G R and C2 > represent some other possible parameters (see examples in 
the Sections 8-13 below), then 

00 

EX,up(0 =C^Y. ^0 = C^iC(l + 2(1 - H)) < 00, (83) 
i=i 

00 

where — ^ Res > 1, is the Riemann zeta-function, and with A^- = A/j, we have 



±2{\t\) 1 

j=l ] = 1 1*1 



where L2 is a slowly varying at infinity function, bounded on every bounded interval. Thus we obtain 
a long range dependence property: 

/ RBup{t)dt = 00. 
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We are going to make an additional assumption that there exists parameters Sj such that 

Ee«^^(") =Ee«^^-^ (85) 

for some random variable Y. The sum 

Y.S,<oo (86) 
must be finite. When we speciahze (81) to this situation we obtain 

oo 

Rsnpit) = C2 5] S, cxp {-A, \t\} , (87) 

for some C2 > 0. This approach allows also to treat the case of several parameters. 
We are starting our considerations with Cq convergence. Let 

corresponds to the process Xj{t). Then, since Xj{-) are independent of each other, 

00 

p(Wl,...,Wg_l) = ]^Pj(Ml,...,Wg_l). (88) 

We have shown above that Pj{ui, . . . , itg_i) is monotone decreasing in ui, . . . , Therefore p{ui, . . . , w^-i) 
being a product of monotone decreasing functions is monotone decreasing itself. Next we prove finite- 
ness of the series. For that recall estimate (67) 



1<P,(6",...,6")< 



VAfo(l)Mo(e-^^^") J - l^Ee^^(o)eE-Y.(o) J 



where C = , and we denote 

Mo(C) = Ee'^'^^^*). 

Then, using (87), we obtain 

1 < p(6",...,6") < C«^r=i'5.'="^^'" ^(j^R...Ab") < i + o(l)-^lnCi?,„„(6"). 

Then X^^i ^(^"7 • ■ • : ^") is finite if the sum 

00 

^i?s„p(6") < 00 

n=l 

is finite. 

We are left to check the mixing condition (20). But this condition follows from the fact that it 
holds for pj, representation (88) and monotonicity of pj. Indeed let 1 < ii < 12 < . . . < *m and put 
^ii,...,i,„ (^) = (wi, . . . , W(j_i), where Ui ~ A if i Cz and otherwise. Let be a number 

which we let to 00 later. Then, for fixed N, 

N N 

i=i j=i 
by the corresponding property of the geometric Ornestein-Uhlenbeck process, see (68). The product, 

00 00 

1 < n PJi^n,...,^AA)) < n pd^'^u...,A^)) ^ 1, 

j=N j=N 
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as A'^ — ^ oo, uniformly in A > 1. Therefore, the mixing property holds. 

Now we turn to proving the scaling property. For that we use Theorem 7. Let A/o(C) = Ee^^ be 
the moment generating function of Y and /(C) = InMo(C). Then, similarly to (69), 



V Mo(g)A/o(e-^.-) 



< 1 - n e-^^^f'^'^^-f^'))'' < .(/'(,) - /(I)) ^a- 

The convergence of the series immediately follows from this estimate and we can apply Theorem 7. 
Note that 

oo 

logEexp{CiX,up(ti) + C2^sup(t2)} = ^logEcxp{CiX,(ti) + C2^j(t2)}, 

where logEexp{Ci-'irj(ii) + C2Xj{t2)},j = 1, 2, .. are given by (57). 
Define the mother process as geometric process 

A(t) = e^-p(*)-=-'^,cx = logEe-^-p(°), A/(C) = Ee^(^-p(°)-'=^\ 

where Xsup(i),t e M, is the infinite superposition process (79). 
Denote 

M(Ci,C2;ii -t2) =exp{-CA-(Ci+C2)}Eexp{CiXeup(ti)+C2^,,p(t2)}. (89) 
We arrive to the following result. 

Theorem 14. Let Xsup(i),i G R+ 6e an infinite superposition of OU-type stationary processes (79) 
such that (80), (82), (85) are satisfied as well as (86). Assume that the Levy measure v in (4-9) of the 
random variable Xsup(i) satisfies the condition that for a positive integer q* € N, 



Then, for any fixed b such that 



xe'^ ^z^(dx) < oo. (90) 

x\>l 



stochastic processes 

l-t n 

Mt) = / n^^'^ {sV)ds,te[Q,l] 
•'° j=0 

converge in Cq to the stochastic process A{t) G Cq,t € [0,1], as n ^ oo. The limiting process A(t) 
satisfies 

EA'?(t) - qe[0,q*]. 

The Renyi function is given by 

T(Q) = g-l-logbEA«(t), [0,(7*], te [0,1]. (92) 

In addition, 

YavA[t)^ [ [ M{Ci,C2;ti-t2)dtds, (93) 
Jo Jo 

where the bivariate moment generating functions M{(^i, (^2]ti — t2) is given by (89). 
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8. Log-tempered stable scenario 

This section introduces a scenario which generahze the log-inverse Gaussian scenario obtained in Anh, 
Leonenko and Shieh (2008). Note that the tempered stable distribution (up to constants) arises in the 
theory of Vershik-Yor subordinator (see Donati-Martin and Yor 2006, and the references therein). This 
section constructs a multifractal process based on the geometric tempered stable OU process. In this 
case, the mother process takes the form A(i) = exp {X (t) — Cx} , where X (t) ,t > is a stationary 
OU type process (51) with tempered stable marginal distribution and cx is a constant depending on 
the parameters of its marginal distribution. This form is needed for the condition EA (t) = 1 to hold. 
The log-tempered stable scenario appeared in Novikov (1994) in a physical setting and under different 
terminology. We provide a rigorous construction for this scenario here. 

We consider the stationary OU process whose marginal distribution is the tempered stable distri- 
bution TS{K,S,'y) (see, for example, BarndorfF-Nielsen and Shephard 2002). This distribution is the 
exponentially tilted version of the positive k— stable law S'(k, S) whose probability density function 
(pdf) is 

= (-1)^-^ sinikn.) £(^2-- (;^)""" ,x > 0. (94) 

fc=i 

The pdf of the tempered stable distribution TS{k, S, 7) is 

C X 1/ K 

TT (x) = 7r{x; n,S,^) = e SK,s{x)e~^'^ , x > 0, (95) 
where the parameters satisfy 

Ate (0,1),(5 > 0,7 > 0. 
It is clear that TS{^,S, 7) = IG {S, 7) , the inverse Gaussian distribution with pdf 

and rS'(i,(5,7) has the pdf 
where here and below 

poo 

K^{z) = / e-^™"''("' cosh{vz)du, z > 0, (96) 



is the modified Bessel function of the third kind of index ly, Rev > 0. 

The cumulant transform of a random variable X with pdf (95) is equal to 

logEe'-"^ = <57-(5(7^-2c)\0<C< (97) 

Note that 

F.X{t) = 2kSj^ , VarA:(i) 4k (1 - k) Sj"^ . 

We will consider a stationary OU type process (52) with marginal distribution TS{k, 5,j). This dis- 
tribution is self-decomposable (and hence infinitely divisible) with the Levy triplet (a, 0, i'), where 

v{du) = b{u)du, 

b{u) ^ 2'^<5— -u-^-^e-^,u > 0. 

r(l - k) 



The BDLP Z{t) in (52) has a Levy triplet (a, 0, i>), with 



v{du) = Xuj{u)du, (98) 

'r(i - k) \u 



w(m)=2«5— ^ - + 2_U-''e-^^,w>0, (99) 
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that is, the BDLP of the distribution TS{k, S, 7) is the sum of a Levy process with density 

r(i- k) 

P{t) 

and a compound Poisson process Z{t) = ^ Zk, which has Levy density 

fe=i 



(7^/2)1-" 
r(l - k) 

meaning tliat the summands Zk have gamma distribution r(l — k, ^-g— )) while the independent Poisson 
process P{t) has rate 6jk. 

B'. Consider a mother process of the form 

A(t) ==exp{X(t)-cx} 

with 



cx 



(57 -(5 (7^ -2) ,7>2'=, (100) 
where X (t) ,t €^ of OU processes with TS{K,dj,^) marginal distributions and covariance function 

Rx (t) = 4k (1 - k) Sj"^ exp {-A \t\} . 
The correlation function of the mother process takes the form 

P(r) = ^^^%i^,7>4^ (101) 
' M(2) - 1 ' ' ' ^ ' 

where 

M(C) = e-^'^^Ee^-^^*), 

and the bivariate moment generating function M(l, 1,t) is given by (59), in which the Levy measure 
D is defined by (98), (99), and cx is given by (100). 
Condition (60) becomes 

1 Ji 

if < < K e (0, 1). Note that M{q) exists if (2q)'" < 7. 
We can formulate the following 

Theorem 15. Suppose that condition B' holds, A > and 

Q = {g:0<g<^,7> max{(2g*)^ 4-^}, n € (0, 1), 5>Q}f^ (0, g*), 
where q* is a fixed integer. Then, for any 

b > exp |-75+ -A_ (^i _ 2g*)' - (7* - 2)''| ,7 > max{(2g*)«, 4^^}, 



the stochastic processes 



pt n 

Mt) = / \{^'^'\sV)ds, i e [0,1] 



converge in Cq to the stochastic process A{t) for each fixed t G [0, 1] as n —> 00 such that, Aq (1) G 
/or g G Q, awrf 

EA''(i)^t^(«)+\ 
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where the Renyi function T (q) is given by 



Moreover, 



VarA{t)^ [ [ [M{l,l:u- w) - l]dudw, 
Jo Jo 

where Af(l,l,r) is given by (57), in which the Levy measure v is defined by (98), (99). 

Theorem 15 follows from the Theorem 11. Note that for k = 1/2 Theorem 15 is an extension of the 
Theorem 4 of Anh, Leonenko and Shich (2008). 

In this partial case we arrive to log-inverse Gaussian scenario where the Renyi function is of the 
form: 

5 



7 - ^rf 



and ^ 

Q = {q:0<g<^,7>2,5>0}n(0,g*) (102) 



if 



I 1-9 1-9 J 



and q* is a fixed integer. 

Note that the set (102) is an extension of the log-inverse Gaussian scenario in the Theorem 4 of 
Anh, Leonenko and Shieh (2008) which is obtained for the set 

Q = {<Z:0<(?< y,7>2,5>0}n[l,2]. 

We can construct log-tempered stable scenarios for a more general class of finite superpositions of 
stationary tempered stable OU-type processes: 

rn 

X„sup(i) = I]^j W,^ e [0,1], 

where Xj(t),j = 1 , . . . , m, are independent stationary processes with marginals Xj (t) ~ TS{k, 5j = 
1, m, and parameters Sj, j = 1, m. Then Xmsup(i) has the marginal distribution TS{k, X^jLi '^i, 7), 
and covariance function 

m 

Rmsup{t) = 4k(1 - k)j~ Z^^j cxp{-Aj \t\} . 

i=i 

It follows from the Theorem 13 that the statement of Theorem 15 can be reformulated for Xmsup{t) 
with 6 = J2]Li , ^ being replaced by \j, and 

M(Cl,C2;tl -^2) = exp{-Cjf(Cl +C2)}Ecxp{ClX^sup(tl) +C2^msup(t2)}, 

where 

m 

log E exp{Cl X™ sup (t 1 ) + C2 sup (t2 ) } = XI ^ ^^P{ ClXjiti)+C2Xj(t2)}, 

and logEexp{CiXj(ti) 4- (2Xj{t2)},j — 1, are given by (57). 

Moreover, one can construct log-tempered stable scenarios for a more general class of infinite 
superpositions of stationary tempered stable OU-type processes: 

00 

X,up{t)=Y.Xj{t),te[0,l], (103) 
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where Xj{t),j = 1, m, ... are independent stationary processes with correlation functions exp {—Xj \t\} , 
and the marginals Xj{t) ~ TS{Hi, Sj, 7), j = 1, m, where parameters Sj = i < < 

l,and Xj = X/j. 

Then Xsup(i) has the marginal distribution TS{k, ^j^ expectation 



and the long-range dependent covariance function 



3 



where L3 is a slowly varying at infinity function, bounded on every bounded interval. It follows from 
the Theorem 14 that the statement of the Theorem 15 remains true with S = J^JLi ^j^ ^-^d 

A/(Ci,C2;ii -t2) =cxp{-cx(Ci+C2)}Eexp{CiX,up(ii)+C2^sup(i2)}, 

where 

cso 

logEexp{CiX,up(ti) + C2^sup(i2)} ^Y.^ogEexp{CiXj{ti) + C2^j(<2)}, 

i=i 

and logEcxp{CiXj(ti) + C2Xj{t2)},j — l,..,m are given by (57) with Levy triplet {a,0~i') given by 
(98) and (99). 

9. Log-normal tempered stable scenario 

This subsection constructs a multifractal process based on the geometric normal tempered stable 
(NTS) OU process. We consider a random variable X = fi + I3Y + \/Y e, where the random variable Y 
follows the TS{k, S, 7) distribution, e has a standard normal distribution, and Y and e are independent. 
We then say that X follows the normal tempered stable law NTS{k, 7, /3, /i, S) (see, for example, 
Barndorff-Nielsen and Shephard 2002). In particular, for k = 1/2 we have that NTS{^,'y, (3, ^, 6) 

is the same as the normal inverse Gaussian law NIG{a, /?, fi, 5) with a = a//?^ + 7^ (see Barndorff- 
Nielsen 1998). Assuming = for simplicity, the pdf of NTS{k, 7, /3, fj., 6) may be written as a mixture 
representation 

00 

p{x; K, 7, 13, 0, S) = -^e^'i+^'' [ z-^^^e-^'-'^+'^''K{z; k, S, j)dz, x G M, (104) 
V 27r J 


where 7r(z; k, S, 7) is the pdf of the TS{k, 6, 7) distribution. We assume that 

e M, 5 >0, 7 > 0, /? > 0, K G (0, 1). 



Let a = \/ PP- + 7-^/" and using the substitution s ~ 1/z we obtain 

00 

p(x;k,7,/3,0,(5) = ^e^'^+'5- / s-^l^e-"^^-"'^^" s)ds, x € R, 

V 27r J 


where Sk..&{x) is given by (94). As an example, the pdf of NTS{^, a, (3, 0, 1) has the form 

3/2 \ 



p(x;i,a,/3,0,l) = (^)ie*>+''^| e-i^^'^^-' U^J 



^/s ds. 
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It was pointed out by Barndorff-Nielsen and Shephard (2002) that NTS{k, 7, /3, /z, S) is self-decomposable. 
Thus, there exists a stationary OU-type process X{t),t > 0, with stationary NTS{K,'y, (3, ^, 6) 
marginal distribution, correlation function 

rx{t)=exp{-X\t\}, 

and 

EX{t) = n + 2Kfi5-f^, 
Vm:X{t) = 2k(57'^ - 4^^^^ (k - 1) . 

We see that the variance can be factorized in a similar manner as in Section 7, and thus superposition 
can be used to create multifractal scenarios with more elaborate dependence structures. 
Also one can obtain the following semi-heavy tail behavior: 

p{x- n, 7, /3, 0, 5) ^ 2-+ifc^7 /^(^ +^) a'^+h |^| — 1 e-"l-l+/'- 

as a; — ±00. 

The cumulant transform of the random variable X{t) with pdf (104) is equal to 

logEef^(*) ^ ^l(: + 5^-5{a' + C)')" , |/3 + CI < « = V/^^+y^. (105) 
The Levy triplet of NTS{k, 7, /3, /i, S) is (a, 0, v), where 

i^{du) = b{u)du, (106) 

bin) = ^a'^+i^^ |^,|-('^+^) K^^.ia \u\)ef-,u& M, (107) 

the modified Bessel function K^{z) being given by (96). 
From (54). (105) and the formulae 

K,{x) = K,i~x), K^,{x) = K,{x), 

—K^{x) ^ A'^(.t) - K^^i{x), 

ax X 

we obtain that the BDLP Z(t) in (46) has a Levy triplet (a, 0, v), with 

i>{du) = Xuj{u)du, 

where 

6 1 k2'^'^^ 
U}{u) = -b{u) - ub'{u) = ^=a'*+2 . 



/2^ r(l-K) 

X {{n ~ i) l^r^'^^^) 1 (a \u\)eP- + l^r^^-^) [_!l±ix,^ (a |«|)e^« 



A-,_ . (a I^De'^-a + A'^+i (« |i.|)e'5"/3]}. (108) 



2" ' ^=+5^ ' ■ ' ' ^ |u 

B". Consider a mother process of the form 

A(t) =cxp{X (t)~cx], 

with 

cx = M + '57 - -5 + 7'/" - (/? + 1)') ^ /? < 

where X [t) ,t ^ [0, 1] is a stationary NTS{k, 7, /3, /i, (5) OU-type process. 
Under condition B", we obtain the following moment generating function 

M (C) = E cxp {C {X (t) -cx)} = e-'''^e^'^+^^-^(^'+^''^~^^+'^^'y , |/3 + C| < a, (109) 
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and bivariate moment generating function 

M (Ci, C2; {h - t2)) = Eexp{Ci - cx) + C2 (^(^2) - cx)} 

= e-^-(«i+«^)Eexp{CiX(ti) + C2 (^(^2)} , (110) 

and Eexp{CiX(ii) + C,2 (-^^(^2)} is given by (57) with Levy measure v having density (108). Thus, 
the correlation function of the mother process takes the form 

M(l,l;r)-1 
P^^^ = M(2)-l ' 

where we assumed that j3 < (7^/'^ — 4)/4. 

Note that as z — > cxj the modified Bessel function of the third kind of index v. 

K.{z) = J\^~"'^-'{^ + + •■•), - > 0, 



Condition (60) now becomes 

/ Me«*" |wr'"+^^ K,+ i {a |M|)e'3"dM < 00, 



if + <a = V/^2+^i/K. 

Theorem 16. Suppose that condition B" holds, A > and 



qeQ = {q:0<q<q* < ^/WTl^- P,I3 < (7'^" " l)/2,/x G M, S>0,Ke (0,1).} , 



where q* is a fixed integer. 
Then, for any 



f , ,5(/32+y/«-(/3 + q*)2)''_5*^(/32+7iA_(/3 + i)2)' 
> exp < -67 H ^ 



l-q* 



the stochastic processes, the stochastic processes 



Mt) = / l[A^'\sV)ds, t e [0,1] 
•^0 J=0 



converge in Cq to the stochastic process A{t) for each fixed t €z [0, 1] as n ^ 00 such that A (1) G Cq 
for q Cz Q, and 

where the Renyi function T (q) is given by 



S 

T{q)= 11-- 



(/32 + 7i/-_(/3 + l)2)« 



log b 



q+ 



Moreover, 

VarA{t)'^ I I [M{l,l;u-w)-l]dudw, 
Jo Jo 

where M is given by (110). 
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Theorem 16 follows from the Theorem 11. 

Note that, for k = 1/2, Theorem 16 is an extension to Theorem 5 of Anh, Leonenko and Shieh 
(2008), which is now extended to present a log- normal inverse Gaussian scenarios with Renyi function: 



5 

Tiq) = I 1 - - 



v//?2+72-(/3 + l)^ 



7 



log b 



q- 



log log b 



where we use the notation NTS{^,j, (3, /i, 5) ~ NIG{a, /3, S, fi), where 7 — \J cP- ~ ^'^ . 

Q={q:0<g<g*< ^fW^ - /3, /3 < (7^ - l)/2, e M, (5 >o} , 



and 



If Xj^(t), j = l,...,m, are independent so that Xj{t) ~ NTS{k,^, (5, fij, Sj), j = l,...,m, then we 
have that 

m m 

+ ... + x,„(t) ~ ivr5(K, 7, /?, E E -^j-)- 

i=i i=i 

We can construct log-normal tempered stable scenarios for a more general class of finite superpositions 
of stationary normal tempered stable OU-type processes: 

rn 

x™,up(t) = 5I^jW'^e [0,1], 

where Xj{t),j = 1 , . . . , m, are independent stationary processes with marginals Xj (t) ~ NTS{k, j, 13, ^j,6j), 
j = 1, m, and parameters Sj, ^J,j,j = 1, ...^m. Then Xmsupit), t G [0, 1] has the marginal distribution 
NTS{k, 7, /3, X^jli Mil X^jli '^j)' ^^'^ covariancc function 

-Rmsup(i) = [2^7"^ - 4k^2 _ ^) j ^ j^. cxp{_Aj |t|} . 

i=i 

EX,„sup(t) = M + 2k/3<57'^, 
VarX„ ,,p(t) = 2k(57'^ - ^k^^S {k - 1) 7^ . 

Then the statement of the Theorem 16 can be reformulated for X^sup with S — X^Jli ^ji ^^'^ 

M (Cl,C2; (tl - t2)) = Eexp{Cl {X„-,sup{tl) - Cx) + C2 {Xmsup{t2) ~ Cx)} 

= e-^-^(';i+«^)Eexp{CiX™,up(ti) + C2 (X,„,up(t2)} , 

where 

ni 

logEexp{CiX™sup(ti) + C2^msup(t2)} = ^logEexp{CiXj(ii) + C2^j(i2)}, 

J=l 

and logEexp{CiXj(ti)-|-(^2-'^j(i2)},J = 1, .., m are given by (57) with Levy measure having density 
(108). 

Moreover, one can construct log-tempered stable scenarios for a more general class of infinite 
superpositions of stationary normal tempered stable OU-type processes: 

00 

Xsup(i) = E^^(*)'*e [0,1], (112) 
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where Xj(t), j = 1, m, are independent stationary processes with correlation functions exp {—Xj \ t\} , 
and the marginals Xj{t) ~ Xj{t) ~ NTS{K,j,l3,fj,j,Sj),j = 1,...,to,,j = 1, m, ... where parame- 
ters S, = = j-iW{i-H))^i <\, and Aj = X/j. 

Then Xsup(i), i G [0, 1] has the marginal distribution NTS{k, 7, /?, X^^i Mji "^i)' expecta- 

tion 

00 00 
i=i j=i 

and the long-range dependent covariance function 

i?.upW = [2^7"^ - 1)7"^] E^,cxp{-A, |t|} = y^'^^) , ^ < ^ < 1' 

where L4 is a slowly varying at infinity function, bounded on every bounded interval. Then the 
statement of Theorem 16 remains true with 5 = X^jli ^jt — X^j^i ^^^^ 

M (Ci, C2; {h -h)) = Ecxp{Ci (X,up(ii) ~cx) + C2 (Xsup(t2) -ex)} 
= e-"^('^^+'^=)Eexp{CiXsup(<i) + C2 (Xsup(t2)} , 

where 

00 

logEexp{CiX,up(ti) + C2^supfe)} = ^logEexp{CiX,(ii) + C2^j(t2)}, 

i=i 

and logEexp{CiXj(ti) + C2Xj{t2)}, j = l,..,m, ... are given by (56) with Levy measure i> having 
density (108). 

10. Log-gamma scenario 

The log-gamma distribution is well-known in the theory of turbulence and multiplicative cascades 
(Saito 1992). In this section, we propose a stationary version of the log-gamma scenario. We will 
use a stationary OU type process (46) with marginal gamma distribution r(/3,a), which is self- 
decomposable, and hence infinitely divisible. The probability density function (pdf) oi X{t),t e M+, 
is given by 

7r(x) = ^x''-ie-"n[o,oo)(a:),a > 0,/3 > 0, (113) 
with the Levy triplet of the form (0, 0, v), where 

u 

while the BDLP Z{t) in (46) is a compound Poisson subordinator 

p(t) 

Zit) = J2 Zn, 

Zn,n = 1,2,..., being independent copies of the random variable r(l,a), and P{t),t > 0, being a 
homogeneous Poisson process with intensity /3. The logarithm of the characteristic function of Z{1) is 



Ac(z)=logEe*^^(i) = ^ 



a — iz 



and the (finite) Levy measure z/ of Z(l) is 

i>{du) = a/3e"""l(o,oo)(w)du. (114) 
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The correlation function is then 



rx (t) ~ cxp (-A \t\ 



B'". Consider a mother process of the form 



k(t) = cxp (X (t) -ex), cx = log ■ 



1 



,a > 1, 



where X (t) ,t g IR+ is a stationary gamma OU type stochastic process with marginal density (113) 
and covariance function 

Rx{t)^^c^p{~X\t\). 



Under condition B', we obtain the following moment generating function: 



M(C)=Eexp(C(X (t)-cx)) = 



1 - 



, C < "7 a > I7 



(115) 



and the bivariate moment generating function is given by the formula (57), in which the measure v is 
given by (114), since 

M (Ci, C2; (ti - ^2)) - Eexp (Ci (X(ti) - cx) + C2 (^(^2) - cx)) 
= e-^-(«i+«^)Eexp(CiX(ii) + C2 (^(i2)) 

/3EjLiOe-"fe'^^l[o,oo)(^j-^) 



-Cx(Cl+C2) 



cxp 



E-=iC,e-^(*^-^)l[o,oo)(t,-s) 



ds 



(116) 



M (Ci, C2; (ti - h)) = exp {-cx (Ci + C2)) 



xexpf / / (exp(M^Oe"^^*^"'^l[o,oo)(ij -s)) -1 



-w(EOe-^^*^-^^l[o.oo)(<, -s))l[-i,il (u)^a/3e-""l(o,oo)(w)dwdsj. (117) 

Thus, the correlation function of the mother process takes the form 

^ M(l,l;r)~l 
M{2) - 1 ' 

where M(2) is given by (115) and M(l, 1; r) is given by (117). It turns out that, in this case. 



log,EA(<)^ 



log 6 



-9 log 



/Slog 1 



and the condition (60) of Theorem 11 holds, since 



ni>i ^ ^ r(/3) 



/-oo 



e« "e-""dM < oo,g* < a. 



We formulate the following 

Theorem 17. Suppose that condition B'" holds and let Q — {q : < q < q* < a,a > 2, P > 0}, 

where q* is a fixed integer. Then, for any 



b > 
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the stochastic processes An (t) defined by (3) for the mother process as in condition B'" converge in 
Cq to the stochastic process A(t) as n ^ oo, such that A (t) € Cq and 

where the Renyi function T (q) is given by 
Moreover, 

VarA(t)^ / / {M{l,l;u-w)-l)dudw, 
Jo Jo 

where M is given by (116) or (117). 

Proof. Theorem 17 follows from the Theorem 11 □ 

Note that the Theorem 17 is an extension of the log-gamma scenario in the Theorem 3 of Anh, 
Leonenko and Shieh (2008) which is obtained for the set 

Q = {(? : < (? < a, a > 2, /3 > 0} n [1, 2]. 

We can construct log-gamma scenarios for a more general class of finite superpositions of stationary 
gamma OU type processes defined in (72), where Xj{t),j = 1, m, are independent OU type gamma 
stationary processes with marginals T{l3j,a),j = l,...,m, and parameters Xj,j = l,...,m. Then 
Xmsup{t),t € M+ has the marginal distribution l^jjCi) S'Hd covariance function 



^ m 



a 

Theorem 17 can be reformulated for the process of superposition X^sup with /3 = X^JLi f^j ^^'^ 
M (Ci, C2; {h - 12)) = Eexp{Ci (X„,up(<i) - cx) + C2 {X^sup{t2) - cx)} 

= e--^(Ci+C2)Ecxp{ClX„,,up(ii) + C2 (Xrnsup(i2)} , 

where 

m 

logEexp{CiXmsiip(ii) -I- C,2X 

m sup (^2)} - ^logEexp{CiX,(ti) + C2^j(t2)}, 
i=i 

and logEexp{Civ'^j(ii) -f C2-''^j"(^2)}, j = 1, •■, rn are given by (116) or (117). 

Moreover, one can construct log-tempered stable scenarios for a more general class of infinite 
superpositions of stationary gamma OU type processes: 

00 

^sup(i) = ^^j(0,te [0,1], 
i=i 

where Xj{t), j = 1, m, .... are independent stationary processes with correlation functions exp {— Aj \t\} . 
and the marginals Xj{t) ~ T{(3j,a),j = 1,...,to,... where parameters /3j = 1 < H < 

l,and Aj = A/j. 

Then Xsup(t), t G [0, 1] has the marginal distribution r(^°^j^ Pj, a), an expectation 



^ 00 

EX,,p(i) = - V /3. 



and the long-range dependent covariance function 

^-pW = ^2 E/^. e-P{-^. 1*1} = ^^T§) . ^ < ^ < 1' 
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where L5 is a slowly varying at infinity function, bounded on every bounded interval. Then the 
statement of the Theorem 17 remains true with /3 = ft ^^'^ 

M {(1X2; {h - t2)) = Ecxp{Cl (X,up(tl) - cx) + C2 (Xsup{t2) - cx)} 

= e-^-'^(^i+«=)Eexp{CiX,,p(ti) + C2 (^sup(i2)} , 

where 

00 

logEcxp{CiX,up(ti) + C2^sup(i2)} = ^logEexp{CiXj(ti) + C2^j(i2)}, 
and logEexp{CiXj(ii) + C2-''^i(^2)}, J = 1, ■■,'m are given by (116) or (117). 
11. Log- variance gamma scenario 

The next example of a hyperbolic OU process is based on the variance-gamma distribution (see, for ex- 
ample, Madan et al. 1998, Finlay and Seneta 2006, Carr et al. 2007). We will use a stationary OU type 
process (46) with marginal variance gamma distribution (A, a, /3, /i), which is self-decomposable, 
and hence infinitely divisible. The pdf of X{t),t £ is 

^(^)- ^ \x-^i^'/'K,^,/2ia\x-^,\)cP^^-^\ xeR, (118) 

VTi-r (K) (2a) 

where K\ (z) is defined by (96) and the set of parameters is 

-f^ = - P^,K> 0,a> >0,neR. 

Note that 

EX(t).,+ ^, VarX(t).= ^(l + 2(^)'), 
and VG {K,a, /3, fi) has semi-heavy tails. The moment generating function of VG {k, a, /3, fi) is given 

by 

logEe«^(*) =fiC + 2k log (^7/\/«2-(/3 + C)') , 1/3 + CI < «■ 

Thus, if Xj (t) j = 1, m, are independent so that Xj ~ VG {Kj, a, /3, fij) , j = I, m, then we have 
that 

Xi{t) + ... + Xmit) ~ VG (ki + ... + Kn, tt, /3, fii + ... + /!„) ■ 

The Levy measure of X{t) has density 

p(u) = Ae^"""'"'," G (119) 
\u\ 

By (54) the Levy measure D of the BDLP Z{t) in (46) has density 

q (u) = (u) - up' (u) , (120) 
_ / -5e"(^+")(/3 + a) + ^e-(^+"), u < 0, 

The correlation function of the stationary process with marginal density (118) is then 

rx{t)^exp {-X\t\). 
B"". Consider a mother process of the form 

A{t) = exp {X{t)-cx),cx=pt + 2k log (^-f/^a^ - (/? + 1)'^ , |/3 + l| <a. 
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where X (t) ,t G R+ is a stationary VG (k, a, (3, ji) OU type process with marginal density (118) and 
covariance function 

i?^ (i) = ^ (^1 + 2 (^^) ^exp(-A|t|). 
Under condition B"", we obtain the moment generating function 

M (C) = E exp (C {X (t) - cx)) = e--xCeK+2«iog(7/V"^-(^+C)^) ^ \/3 + C\ < a, (121) 

and the bivariate moment generating function 

M (Ci, C2; (ti - t2)) = E exp (Ci (X(ti) - cx) + C2 iXit2) - cx)) 

= e-^-'^e^i+^^^Ecxp (CiX(ti) + C2 (^(^2)) , (122) 

where Eexp {(iX{ti) + C2 (-''^(^2)) is given by (57) with Levy measure having density (120). Thus, 
the correlation function of the mother process takes the form 

M(l,l;r)-1 
Af (2) - 1 ' 

where M(2) is given by (121) and A/(l, 1; r) is given by (122). 
The condition (60) of Theorem 11 holds for g < a — |/3| . 

Theorem 18. Suppose that condition B"" holds and let 

Q = {q:0<q<q* <\a\-\p\,\p + l\<a,K>0}, 

where q* is a fixed integer. 
Then, for any 



b > exp < 2k 



1 7 

1-9* ^ ^a^-{l3 + q*)^ 



1 1 
■log 



q* v/a'-(/3 + l)' 



the stochastic processes A„ [t) defined by (3) for the mother process as in condition B"" converge in 
Lq to the stochastic process A(t) as n 00 such that, if A (1) G Lq for q G Q, 

where the Renyi function is given by 

^(g)=g(i + T:^iog^^ T. ...J + log \/"' - (/3 + 9)' - log 7 - 1- 



Moreover, 



log6 - (/? + 1)2 y log6 iQg^ 

VarA(t)^ / / {M{l,l;u-w)-l)dudw, 
Jo Jo 



where M is given by (122). 

Proof. Theorem 18 follows from the Theorem 11. □ 

We can construct log- variance gamma scenarios for a more general class of finite superpositions 
of stationary variance gamma OU type processes of the form (72), where Xj{t),j = l,...,m, are 
independent variance gamma stationary processes with marginals 

Xj ~ VG {Kj,a, P, fij) ,j = 1, ...,m, and parameters Xj,j = 1, ...,m. Then Xmsup(t),t G M+ has the 
marginal distribution VG {ki + ... + Km,a, f3,S, fii + ... + /.(m) and covariance function 

Rmsnpit) = ^ 1^1 + 2 (^^^ ^^K,exp(-Aj |t|),teR. 
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The generahzation of Theorem 18 remains true for this situation with A = SJLi ^ fJ- ~ SJLi 
and 

M {C1X2; ih - t2)) = Ecxp{Cl (X^snpih) - Cx) + C2 (^msup(i2) " Cx)} 

= e-=^(^i+«^)Eexp{CiX™,up(ti) + C2 (^msup(t2)} , 

where 

m 

logEexp{CiX™,up(<i) + C2^msup(t2)} = ^logEexp{CiXj(ii) + C2^j(i2)}, 

and logEexp{CiXj(ti) + C,2X , j — 1, --^fn are given by (122). 

Moreover, one can construct log- variance gamma scenarios for a more general class of infinite 
superpositions of stationary gamma OU type processes: 

00 

^sup(i) G [0,1], 

where Xj (t), j = 1, m, ... are independent stationary processes with correlation functions exp {— A^ \t\} , 
and the marginals (t) ~ VG {Kj , a, fij) , j = 1, m, ... where parameters kj = i^ij = 

i <H< l,and A, = A/j. 

Then Xsupit),t G [0, 1] has the marginal distribution VG ^X^jli o^, Mj) , ^-ii expectation 

and the long-range dependent covariancc function 

Rs.At) - :^ (1 + 2 ± cxp{-A, |t|} ^g^, i < < 1, 

where Lg is a slowly varying at infinity function, bounded on every bounded interval. Then the 
statement of the Theorem 15 remains true with fi ~ I^JLi Mj, = Sj=i ^"^^ 

M (Ci, C2; (ti - is)) = E exp {Ci (X,,p(ti) - cx) + C2 (^sup(t2) ~ cx)} 
= e-^^(^i+«=)Eexp{CiX,,p(ii) + C2 (^sup(i2)} , 

where 

00 

log E cxp{CiX,up (ti ) + C2^sup (^2 )} = ^ log E exp{Ci X^ {h ) + C2X, (is ) } , 
and logEexp{CiXj(ii) + (2Xj{t2)}J = 1, .., ... are given by (122). 
12. Log-Euler's gamma multifractal scenario 

This section presents a new scenario which is based on Euler's gamma distribution (see, for example, 
Grigelionis 2003). We consider the random variable Y with the gamma distribution r(/3, a) having 
pdf (113) and the random variable 

X^=7logr,7 7^0, 

which has the pdf 

where the parameters satisfy 

a > 0,/3 > 0,7 ^ 0. 



D. Denisov and N. Leonenko / Limit Theorems for Multifractal Products 



39 



The characteristic function of random variable X with pdf (123) is 

r(/3)a'T^ 

Grigehonis (2003) proved that for 

J >0,a>0,^>0,77^0, (124) 

the function 



r(/3 + i-fz) 



s 



is a self-decomposable characteristic function. We denote the distribution of the random variable X 
byr(7,«,/3,5). 

We note that r(7, , 1, 1), 6* e R, is the Gumbel distribution with location parameter 6 and scale 
parameter I7I , since with 

A(a;) = exp {-e"""} , A{x) = 1 - A(-a-), xER, 

( a(^), 7<0, 
P{X<x}={ ,x€R. 
\ a(^), 7>0 

We will use a stationary OU-type process (46) with marginal distribution r(7,a,/3, 5), which is self- 
decomposable, and hence infinitely divisible. It means that the characteristic function of X{t),t G 
is of the form (125) under the set of parameters (124). Note that, for /3 > 0, we have 

r(/3 + iz) = r(/3) expjiz J ^ l{o<x<i}^ dx 





+ J (e"" - 1 - zzxl{_i<,<o}) _ ^,^ dxj, 

— oc 

and thus the distribution corresponding to the characteristic function (125) has the Levy triplet 
((5a, 0, v), where 



and 



a = -f I (— _^ _ 1 + 7 / ^ — dx-7loga, 



l7l 



v{du) = Sb{u)du, (126) 
^T-Y, u<0,7>0, 



1, 

e -I 



bin) 



\u\(l 



, it>0,7<0. 

l-ei" 



Thus, if Xj{t), j = 1, m, are independent so that Xj{t) ~ r(7, a, /3, Sj), j = 1, m, then we have 
that 

Xi{t) + ... +X„(t) ~ r(7,a,/3,5i + ... + Sm) 
and if j = 1, m, are independent so that Xj{t) ^ r(7, a^, ^, (5), .7 = 1, m, then 

Xi(t) + ... + X„(f) ~ r(7, n a„ /3, 5). 
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The BDLP Z(t) in (52) has a Levy triplet (a, 0, i>), where 

.ir(/3) 



40 



r(/3) 



+ jXS\oga — XS / xLo{x)dx, 



a: >1 



with the density of v given by 



v[du) = \5u!{u)du, 



(127) 



uj{u) = < 



I 1 - e 



— e T 1 — e 



|7 



7 > 0,u < 0, 



V-r?, 7<0,w>0. 



l-e 



The correlation function of the stationary process X (t) then takes the form 

rx (f) =exp{-A|t|}. 

Note that 



— r (B) 

EX(t) = 7^^Y^ " 7'^loga, VarX(t) - ^7^ 



-fix 



1 - e^^ 



-da;. 



B^. Consider a mother process of the form 



with 



A{t) = exp {X (t) ~cx}, t e R+ 



rf/3 + 7) 

cjc=^log ^7^, y ,/3>0,7<0,/3>-7. 



r (/3) e7 

where X (t) , t G [0, 1] is a stationary r(7, a, (3, 6) OU-type stochastic process with covariance function 

Rx (t) = YarX (t) exp {-A . 
The logarithm of the moment generating function of r(7, a, /3, S) is 

logEe'^^^*) = Slog^^^,0 < C < > 0,7 < 0. 

r (/3) eTC 7 

Under condition B^, we obtain the following moment generating function 
M (C) = Eexp {C {X (t) - cx)} = e-^-^e''^<^ , < C < 

7 

and bivariate moment generating function 

M {C1X2; ih - t2)) = Eexp{Ci (Xih) - cx) + C2 iX{h) - cx)} 



(128) 



_ p-Cx(Cl+C2) 



Eexp{CiX{h) + C2{X{t2)}., 



(129) 



where Eexp{CiX(ti) + C2 (-''^(^2)} is given by (57) with Levy measure i> having density (127). Thus, 
the correlation function of the mother process takes the form 



Pir) 



A'/e(l,l;r)-l 



Me (2) - 1 ' 

where Mo{2) is given by (128) and Me(l, 1; r) is given by (129). 



(130) 
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Theorem 19. Suppose that condition holds and 

qeQ=\q:0<q<q*< --,13 > 0, 7 < 0, ^ > -7, (5 > 0, a > 
I 7 

where q* is a fixed integer. Then, for any 

"r*(/3 + 7'z*)r«*(/3) ~ 



41 



b > 



r(/3)r9* (/3 + 7) 



the stochastic processes 



An{t)^ / TTA^J") e [0,1] 



converge in Cq to the stochastic process A(t), t e [0, 1] as n — > 00 such that A{1) € Cq, and q Q, 

EA'^{t)^e^''^+\ 

where the Renyi function is given by 

6 



T(g)=,(l + ^logr(/3 + 7)-^logr(/3) 



log 6 



logr(/3 + g7)-l,geQ. 



Moreover, 



VarA{t)^ / / [Mg{l,l;u- w) - l]dudw, 
Jo Jo 



where Mg is given by (129). 

Proof. Theorem 19 follows from the Theorem 11 



□ 



We can construct log-r(7, a, /3, (5) scenario for a more general class of finite superpositions of sta- 
tionary OU-type processes: 

m 

^msup(i) = II^jW,ie [0,1], 
where Xj{t),j = 1, ...,m, are independent stationary processes with marginals Xj{t) ~ r(7, a, /?, 

m 

j = 1, m, and parameters Sj, j = 1, m. Then Xmsup{t) has the marginal distribution r(7, a, /3, ^ (5 
covariance function 



3/' 



Rsup{t) 



7 / X 



1 - e- 



-dx 



^J,cxp{-A,|i|} 



and 



M (Ci,C2; {ti - t2)) = Eexp{Ci(X„,up(ii) - cx) + C2(X„,sup{t2) ~ cx)} 
= e-^^(^i+f^)Eexp{CiX™,,p(ii) + C2 (^msupfe)} , 



where 



logEcxp{CiXmsiip(ti) + C2X 

m sup (^2)} = ^logEcxp{CiX,(ti) +C2X,(t2)}, 

and logEcxp{(iXj(ii) + C2Xj{t2)},j — 1, .., m arc given by (129). 

Moreover, one can construct log-r(7, a, /3, (5) scenarios for a more general class of infinite superpo- 
sitions of stationary gamma OU type processes: 



(131) 
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where Xj (t), j = 1, m, ... are independent stationary processes with correlation functions exp {—Xj \t\} , 
and the marginals Xj{t) ~ T{'j,a, (3,Sj),j = 1, ...,m, ... where parameters Sj = ^ < H < 

l,and Xj = X/j. 

Then Xsup{t),t G [0, 1] has the marginal distribution r(7, a, /3, X^j^i ^j)^ varianccVarX(t) ~ Sj^ 
expectation 

— r(/3) °° 

EX{t) = {j^^^-j\oga)J2S„ 
and the long-range dependent covariance function 

^sup(t) = 7^ / ^^J-Y.^, exp {-A, \t\] = 2 < < 1, 

"'0 ' j=l I'l 

where L7 is a slowly varying at infinity function, bounded on every bounded interval. Then the 
statement of the Theorem 19 remains true with S = ^^'^ 

M (Ci, C2; ih - ia)) = Eexp {Ci (X.^pih) - cx) + C2 (^sup(t2) - cx)} 
= e-^^(^i+«^)Eexp{CiX,,p(<i) + C2 (Xsup(t2)} , 

where 

00 

log E exp{Ci X,up (i 1 ) + C2 Xsup (<2 ) } = ^ log E exp{Ci X^ (t 1 ) + C2^j (^2 ) } , 
and logEexp{CiXj (ti) + C2^j (t2)}, j = 1, ■•, •■• are given by (129). 



13. Log-z scenario 

The next scenario is based on the z-distribution (see, for example, Grigelionis 2001). We consider a 
pdf of the form 

27rexp(2^(x-/.)) 

7r(a;) = ^ ^ ^, xeR, (132) 

aB(/3i,/32)(l + exp(2^(.-M)))'^'''^ 

where the set of parameters is 

a > 0,/3i > 0,^2 > 0,/i e M 

(see Prentice 1975, Barndorf-Nielsen et al. 1982). The characteristic function of a random variable X 
with pdf (132) is given by 

T^jzx _ Bjfii + ^,/32 - It) , c 

- — w^) — ' 

This distribution has semi-heavy tails and is known to be self-decomposable (Barndorf-Nielsen et al. 
1982), hence is infinitely divisible. Due to this infinite divisibility of the z-distribution, the following 
generalization can be suggested. 

We will use a stationary OU type process (46) with marginal generalized z-distribution Z{a, /32, 6, fi). 
The characteristic function oi X{t), t £ R+ is then of the form 

C iaz iaz \ 

— — ) ' ' ^ ^ 

where the set of parameters is 

a > 0, ^1 > 0, ^2 > 0, (5 > 0, A* € K. 
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This distribution is self-decomposable, hence infinitely divisible, with the Levy triplet (a,0, z/), where 



and 



1 - e-^ 
v{du) = b{u)du, 

2Se ~^ 



-dx + ^, 



(134) 



b{u) 



1— e a 



2Se' 



lu\[l-e 



u> 0, 



u<0. 



Thus, ii Xj{t), j = 1, m, are independent so that Xj{t) ~ Z{a, /3i, /32,Sj, fij), j = 1, m, then we 
have that 

The BDLP Z{t) has a Levy triplet (a,0, where 



a = Xfi + 



dx — \ xuj{x)dx, 



\x\>l 



with the density of being given from 



iy{du) = \u){u)du, 



^ (/32e- 



1 - e- — " ) +e 

2T(f)i+l) 



2^(32 + 1) . 



x2 , W > 0, 



u < 0. 



(135) 



The correlation function of the stationary process with marginal density (132) is then 

rx (t) =cxp(-A|t|). 
The pdf of the generalized z-distribution Z{a, /32, (5, fi) has semi- heavy tails: 

7r(x) ~ C± Ixl"* e-°'±l^l, |x| ^ ±oo, 

where 



p± = 2d - 1,(7+ = = ,C± 

Q' a 



27r 



Note that 



EX{t) = 



1 - e-^' 



dx + /I, VarX(i) = 



aB(/3i,/32)y r(2,5) 



(27r)2 y„ 1-6-=^ 



In particular, the generalized z-distribution Z{a., 5 + ^15 — ^:^>/^) = Af (a, /3, (5, \i) is known as the 
Meixner distribution (Schoutens and Tcugcls 1998, Grigelionis 1999, Morales and Schoutens 2003). 
The density function of a Meixner distribution is given by 



TT [x] 



(2cosf) 



2(5 



cxp ( — (x — fl) 



T{d + i 



X — fl 



where 
Note that 



27rar(2d) ^ V" 

a > 0, -TT < /? < TT, (5 > 0, At e M. 
\T{x + ty)\^ ~ ^/2^|y|"="^/2e-'^l''l/2 as \y\ 00. 



, X e 
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This distribution is infinitely divisible and self-decomposable with triplet (a, 0, v) , where 

P p sinh(/?x/2) ^ 

a = ao tan — — Id I . , , rrrax + M 

2 sinh(7ra:/2) ^ 

and 

<5expf^ 

The cumulant function is 



V (du) — -, — -T-du. 

^ ' usinhf^) 



cosh ((az — ip) /2) 

In particular, the hyperbolic cosine distribution Z{a, i, i, 0, ^) — M {a, 0, i, /x) has the pdf 

, , 1 

a cosh (-(x - 

and characteristic function 

cosh (if)' 

while the logistic distribution Z{a, 1,1,0, ^) has the pdf 

27rexp (-(x — /i)) 
a(l + cosh(f(.T-M))) 

and characteristic function 

g »^X(t) ^ ^^^^ g 

2sinh(^) 

Another example is the z-distribution Z(27r, 0, log |^), which is the log -Ffci.fcj distribution, 

where Fk^^k2 is the Fisher distribution (Barndorff-Nielsen et al. 1982). Note that the generalized z- 
distributions and generalized hyperbolic distributions form non-intersecting sets. However, one can 
show that some Meixner distributions and corresponding Levy processes can be obtained by subor- 
dination, that is, by random time change in the Brownian motion (see, for instance. Morales and 
Schoutcns 2003). 

B^^. Consider a mother process of the form 

A(t) =exp(X it)^cx), 

with 

c.^2s( logr + ^)+ logr -^)- log ^ (^^^ 



where X (t) , t £ is a stationary Z(a, /3i, /32, S, fi) OU-type process with covariance function 

Rx (t) = (VarX (t)) cxp (-A \t\) . 
The logarithm of the moment generating function of Z{a, /32, S, /i) is 

iogEeC-(*) = 2s(io,r(p, + ^)+ logr log + K, Ce^ '^^^ '^^^ 



Under condition B^/, wc obtain the moment generating function 



M(C) =Eexp(C(X(t)-cx)) = e-^^^e^^^'^(*», (136) 
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and the bivariate moment generating function 

M (Ci, C2; [h - t2)) = E exp (Ci {X{h) - cx) + C2 iX{t2) - cx)) 

= e-«(Ci+C2)Ecxp (CiX(ti) + C2 (^(^2)) , (137) 

where Eexp (CiX(ti) + C2 (-^^(^2)) is given by (56) with Levy measure i> having density (135). Thus, 
the correlation function of the mother process takes the form 

M(l,l;r)-1 
- Af(2)^l ' 

where M (2) is given by (136) and M(l, 1; r) is given by (137). 

Theorem 20. Suppose that condition B^^ holds and let Q = |g e (0, q*) : < ^ < q* < /3i < 13i + ^ > 0, ^2 
where q* is a fixed integer. 
Then, for any 



b> 



-, 25 

r(/3i)/r(/?2) 
r(A + ^)r(/32-^) 



the stochastic processes An (t) defined by (S) for the mother process (95) converge in Cq to the stochas- 
tic process A{t) as n —> 00 such that A (1) G Cq for q Q, 

EA{ty ^ t^^''^+\ 

where the Renyi function is given by 



T{q)^q \ 1 + 



26 (logr (/3i + ^) + logr {P2 - If) - log — 



Moreover, 



log 6 

VarA(t)^ / / {M{l,l;u-w)-l)dudw, 
Jo Jo 



where M is given by (137). 

Proof. Theorem 20 follows from the Theorem 11. □ 

We can construct log-z scenarios for a more general class of finite superpositions of stationary OU- 
type processes of the form (72), where Xj(t),j = 1, ...,m, are independent stationary processes with 
marginals 

Xj(t) ~ Z{a, l3i, l32,Sj, ^ij), j = l,...,m and parameters Sj,^j, j = l,...,m. Then Xmsup{t),t G ^+ 

m m 

has the marginal distribution Z{a, Pi, P2, S ^ji S f-j) covariance function 

/ 2^2 f°° e-^'^^ + e~l^^^ \ -A 
i?msup(i) = \J^^ ^ 1 _ ^-x ^^j 2^^SjC^vi->^3 \A) ■ 

The statement of Theorem 20 can be reformulated for the process of superposition X^ sup {t) with 
5 Sj , fi^ Mj, and 

M (Cl, C2; iti-t2))^ Eexp{Cl (X„sup(tl) -Cx)+C2 (^r„sup(i2) - cx)} 
= e-^^(^i+«^)Eexp{ClX,„sup(tl) + C2 (^rnsup(i2)} , 
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where 

m 

logEexp{ClX„,up(ii)+C2^msup(t2)} =X!l°gE^^P{Cl^j(il) +C2^j(i2)}, 

J = l 

and \ogEexp{(iXj{ti) + C2Xj{t2)},j = 1, ..,m are given by (137). 

Moreover, one can construct log-Z scenarios for a more general class of infinite superpositions of 
stationary OU type processes with marginal z-distribution: 

oo 

^.upW = e [0,1], 

where Xj{t),j = 1, m, ... are independent stationary processes with correlation functions exp {—Xj |t|} , 
and the marginals ~ Z{a, f3i, (32,Sj, — 1, to, ... where parameters (5j ^ i < 

H < l,and Xj = X/j. 

oo oo 

Then Xsup(t),t e [0, 1] has the marginal distribution Z(q!, /3i, /32, '^j, expectation 



EX,„p(t) = (- / — dx) 2^Sj + }_^ fij, 



and the long-range dependent covariance function 



Rs^At) - ((^ / ^ c^.)^^,-exp{-A, K|} ^1^, 2 < ^ < ^' 



where Lg, is a slowly varying at infinity function, bounded on every bounded interval. Then the 
statement of Theorem 20 remains true with S — J2'jLi ^j-l^ = Xjli Mi ^^d 

M (Cl,C2; ih - t2)) - Eexp{Cl (X,up(tl) - cx) + C2 (^sup(t2) - cx)} 

= e-^-'^(^i+';^)Eexp{CiX,„p(ti) + C2 (^sup(i2)} , 

where 

oo 

logEexp{CiX,up(ti) + C2X,up(t2)} = ^logEexp{CiX,(ti) + C2X,{t2)}, 

and logEexp{CiXj(ti) + C2^j (i2)}, J 1, ■•, •■• are given by (137). 

In principle, it is possible to obtain log-hyperbolic scenarios for which there exist exact forms of 
Levy measures of the OU process and the BDLP Levy process; however some analytical work is still 
to be carried out. This will be done elsewhere. 
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